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We will review our understanding of non-abelian gauge theories in finite physical 
volumes. It allows one in a reliable way to trace some of the non-perturbative 
dynamics. The role of gauge fixing ambiguities related to large field fiuctuations 
is an important lesson that can be learned. The hamiltonian formalism is the 
main tool, partly because semiclassical techniques are simply inadequate once the 
coupling becomes strong. Using periodic boundary conditions, continuum results 
can be compared to those on the lattice. Results in a spherical finite volume will 
be discussed as well. 



1 Introduction 

We have decided to take this opportunity of contributing to a handbook of 
QCD, in honor of Boris Lazarevich loffe, to bring together resuhs and methods 
that were developed in solving for the low-lying spectrum in a finite volume. 
The emphasis will be on the dynamical aspects of the classically scale invariant 
theory for non-abelian gauge theories in 3-1-1 dimensionslil The challenge is to 
understand how the mass gap is generated. Due to the need for regularization, 
at the quantum level breaking the scale invariance, a running coupling appears. 
The difficulty lies in the fact that this coupling increases at low-energies and 
long-distance scales, beyond the point where one has control over the field 
fluctuations, which become so large that they probe the essential non-linearities 
of the theory. A finite volume explicitly breaks the scale invariance. However, it 
does so in a rather mild way. Classically one can use the scale transformation to 
go from one physical volume to another, and only the running of the coupling 
constant prevents us from taking the infinite volume limit. As the longest 
distance scale, i.e. the lowest energy scale, is set by the volume, one can keep 
in check the growth of the running coupling constant. 

We will describe the results mainly in the context of a hamiltonian pictureH 
with wave functionals on configuration space. Although rather cumbersome 
from a perturbative point of view, where the covariant path integral approach 
of Feynman is vastly superior, it provides more intuition on how to deal with 
non-perturbative contributions in situations where semi-classical techniques 
can no longer be used, like for observables that do not vanish in perturbation 
theory. The high energy modes can be well-approximated by a harmonic oscil- 
lator contribution to the wave functional. In the direction of these field modes 
the potential energy rises steeply. Their contribution, including regulating the 
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ultraviolet behavior, is treated perturbatively, giving in particular rise to the 
running of the coupling. The finite volume allows us to have a well-defined 
mode expansion. Due to the classical scale invariance, the hamiltonian can be 
formulated in terms of dimensionless fields. This can be extended to the quan- 
tum theory, as Ward identities allow for a field definition without anomalous 
scaling. Apart from the overall scaling dimension of the hamiltonian, only the 
running coupling introduces a non-trivial volume dependence. 

Due to the non-abelian nature of the theory the physical configuration 
space, formed by the set of gauge orbits A/Q {A is the coUectioiuof connec- 
tions, Q the group of local gauge transformations) is non-trivialQ Most fre- 
quently, coordinates of this orbit space are chosen by picking a representative 
gauge field on the orbit in a smooth and preferably unique way. Linear gauge 
conditions like the X/andau or Coulomb gauge suffer from what is known as 
Gribov ambiguities □ The region of field space that contains no further gauge 
copies is called a fundamental domain for non-abelian gauge theoriesu 

Having arranged the low-energy field modes (and all those modes not af- 
fected by the cutoff) to be scale invariant, the spreading of the wave functional 
is completely caused by an increasing coupling. This is what leads to non- 
perturbative effects. Asymptotic freedom on the other hand guarantees that 
in small volumes the running coupling is small and it thus keeps the wave func- 
tional localized near the classical vacuum manifold. In a periodic geometry, 
this pertjirbative analysis was pioneered by BjorkenQ and further developed by 
LiischerQ. The essential ingredient we have added to address non-perturbative 
effects is boundary conditions in field space, at the boundary of the fundamen- 
tal domain, with gauge invariance implemented properly at all stages. 

The non-perturbative features due to spreading of the wave functional can 
be followed out to a physical volume of about one cubic fermi (setting the 
scale by the string tension or lowest glueball state). This is, on the basis of a 
comparison with lattice Monte Carlo data, the point where in the pure gauge 
theory the confining string is being formed, with no significant finite volume 
dependence beyond this volume. What has become clear is that the transition 
from the finite to the infinite volume is driven by field fluctuations that cross 
the barrier which is associated with tunneling between different classical vacua. 
This is natural, since this barrier (the finite volume sphaleron), will be the 
direction beyond which the wave functional can first spread most significantly, 
as it provides the lowest mountain pass in the energy landscape. 

In Sec. H we describe the process of complete gauge fixing, and the fact 
that the boundary of the fundamental domain, unlike its interior, has gauge 
copies that implement the non-trivial topology of configuration space. This is 
applied in Sec. to formulating non-abelian gauge theories on a torus, both 
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for SU(2) and in Sec. 3.1 for SU(3). The influence of massless quarks on the 
smafl volume vacuum structure is discussed in Sec. whereas Sec. 3.3 gives 
a short review of the non-perturbative evaluation of the running renormal- 
ized coupling, making explicit use of the finite volurae geometrycl In Sec. 3.4 



we consider the case of twisted boundary conditionsp and in Sec. 3.5 we dis- 
cuss supersymmetric Yang-Mills theory in a finite volume in the light of some 
recent developments concerning the Witten indexEJ Sec. || outlines what is 
known about instantons on the torus, reviewing the Nahm transformatiorcJ 
in Sec. 4.1, Sec. ^ analyzes the situation for a spherical geometry, particularly 
well adapted to study the role of instantons in the low-lying glueball spectrum. 
Finally, we shall consider the behavior in large finite volumes. We discuss how 
the previous results fit together and what this implies. In Sees. |6.l| and 



we review the volume dependence whenever the polarization clouds are well- 
contained in the finite volumetHI In Sec. 6.3 we briefiy mention the situation 
in the absence of a mass gap, particularly relevant for QCD wi^h its light up 
and down quarks, for which chiral perturbation theory appliesEj_We conclude 
with a review of 't Hooft's electric-magnetic duality on the toruaa in Sec. 6.4. 



2 Complete Gauge Fixing 

An (almost) unique representative of the gauge orbit is fflwd by minimizing 
the norm of the vector potential along the gauge orbitiJ'Eil 

GA{h) = \\[h]Af = - f tr ((/i-i(x)A,(x)/i(x) +/^-i(x)a,/^(x))') , (1) 

where the vector potential is taken anti-hermitian, the integral over the finite 
spatial volume M is with the appropriate canonical volume form and ft.(x) is 
a Lie-group element, with a short-hand notation for the associated gauge 
transformation. Note that in these conventions the field strength is given by 

F^,(x, t) ^ a^A,(x, t) - 9,A^(x, i) + [A^(x, t), A,(x, i)\ (2) 

and the action by 

f 5"'tr(^^^i(x,i)). (3) 

JM 

Expanding around the minimum of Eq. (|l|), writing h{x) = exp(X(x)) {X{x) 
is, like the gauge field Ai{x), an element of the Lie-algebra) one easily finds 

\\[h]Af = \\A\f + 2[ tT{Xd^A,)+ [ tr {X''FP{A)X) (4) 
Jm Jm 

+ 1 f tY {X[[A,,X],d,X]) + ^ f tr {[V,X,X][d,X,X]) + 0{X'), 
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where FP{A) is the Faddeev-Popov operator {ad{A)X = [A,X]) 

FP{A) = ~d{D,{A) = -9,(3, + ad(A,)). (5) 

At a local minimum the vector potential is therefore transverse, diAi — 0, 
and FP{A) is a positive operator. The set of all these vector potentials is by 
definition the Gribov region fi. Using the fact that FP{A) is linear in A, 
is seen to be a convex subspace of the set of transverse connections F. Its 
boundary dO. is called the Gribov horizon. At the Gribov horizon, the lowest 
non-trivial eigenvalue of the Faddeev-Popov operator vanishes, and points on 
dVi are associated with coordinate singularities. Any point on can be seen 
to have a finite distance to the_nxigin of field space and in some cases even 
uniform bounds can be derivedtSEJ 

The Gribov region is the set of ^oca/ minima of the norm functional, Eq. 
and needs to be further restricted to the absolute minima to form a fundamental 
domain,El which will be denoted by A. The fundamental domain is clearly 
contained within the Gribov region. To show that also A is convex, note that 

II [h]Af -\\Af^ J tr (A^) - y tr + h-'d,hf^ (6) 

= y tr {h^FPf{A)h) = {h,FPfiA)h) , FPf{A) = ^d,D,{A), 

where FPf{A) acts on the fundamental representation and is similar to the 
Faddeev-Popov operator. Both FP{A) and FPf{A) are hermitian operators 
when v4 is a critical points of the norm functional, i.e. for A transverse. We 
can define A in terms of the absolute minima over h E Q oi {h, FPf{A) h) 

A = {A e F|min(/i,FP/(A)/i) = 0}. (7) 

Using that FPf{A) is linear in A, the convexity of A is automatic: A line 
connecting two points in A lies within A. 

If we would not specify anything further, as a convex space is contractible, 
the fundamental region could never reproduce the non-trivial topology of the 
configuration space. This means that A should have a boundaryliJ Indeed, as 
A is contained in fJ, this means A is also bounded in each direction. Clearly 
A = is in the interior of A, which allows us to consider a ray extending out 
from the origin into a given direction, where it will have to cross the boundary 
of A and fi. For any point along this ray in A, the norm functional is at its 
absolute minimum as a function of the gauge orbit. However, for points in U, 
that are not also in A, the norm functional is necessarily at a relative minimum. 



4 



The absolute minimum for this orbit is an element of A, but in general not 
along the ray. Continuity therefore tells us that at some point along the ray, 
this absolute minimum has to pass the local minimum. At the point they 
are exactly degenerate, there are two gauge equivalent vector potentials with 
the same norm, both at the absolute minimum. As in the interior the norm 
functional has a unique minimum, again by continuity, these two degenerate 
configurations have to both lie on the boundary of A. 

When L denotes the linear size of M, we may express the gauge fields 
in the dimensionless combination of LA (in our conventions the fields have 
no anomalous scale dependence), and the shape and geometry of the Gribov 
and fundamental regions are scale independent. We should note that the norm 
functional is degenerate along the constant gauge transformations and indeed 
the Coulomb gauge does not fix these gauge degrees of freedom. We simply 
demand that the wave functional is in the singlet representation under the 
constant gauge transformations and we have A/G = A/Q. Here A is assumed 
to include the non-trivial boundary identifications that restore the non-trivial 
topology oi A/Q. 
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Figure 1: Sketch of the fundamental (shaded) and Gribov regions. The dotted Unes indicate 
the boundary identifications. 

If a degeneracy at the boundary is continuous, other than by constant 
gauge transformations, one necessarily has at least one non-trivial zero eigen- 
value for FP(A) and the Gribov horizon will touch the boundary of the fun- 
damental domain at these so-called singular boundary points. We sketch the 
general situation in Fig. ^. In principle, by choosing a different gauge fixing 
in the neighborhood of these points one could resolve the singularity. If sin- 
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gular boundary points would not exist, all that would have been required is 
to complement the hamiltonian in the Coulomb gauge with the appropriate 
boundary conditions in field space. Since the boundary identifications are by 
gauge transformations the boundary condition on the wave functionals is sim- 
ply that they are identical under the boundary identifications, possibly up to 
a phase in case the gauge transformation is homotopically non-trivial. 

Unfojiunately, one can argue that singular boundary points are to be 
expectcdIiZI Generically, at singular boundary points the norm functional un- 
dergoes a bifurcation moving from inside to outside the fundamental (and 
Gribov) region. The absolute minimum turns into a saddle point and two lo- 
cal minima appear, as indicated in Fig. ^j. These are necessarily gauge copies 
of each other. The gauge transformation is homotopically trivial as it reduces 
to the identity at the bifurcation point, evolving continuously from there on. 



Figure 2: Sketch of a singular boundary point due to a bifurcation of the norm functional. 
It can be used to show that there are homotopically trivial gauge copies inside the Gribov 
horizon (H). 

Also Gribov's original arguments for the existence of gauge copiesi (show- 
ing that points just outside the horizon are gauge copies of points just inside) 
can be easily understood from the perspective of bifurcations in the norm 
functional. It describes the generic case where the zero-mode of the Faddeev- 
Popov operator arises because of the coalescence of a local minimum with 
a saddle point with only one unstable direction. At the Gribov horizon the 
norm functional locally behaves in that case as X^, with X the relevant zero 
eigenfunction of the Faddeev-Popov operator. The situation sketched in Fig. || 
corresponds to the case where the leading behavior is like X^. 

The necessity to restrict to the fundamental domain, a subset of the trans- 
verse gauge fields, introduces a non-local procedure in configuration space. 
This cannot be avoided since it refiects the non-trivial topology of this space. 
We stress again that its topology and geometry is scale independent. Ho- 
motopical non-trivial gauge transformations are in one to one correspondence 




(a) 



(b) 
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with non-contractible loops in configuration space, which give rise to conserved 
quantum numbers. The quantum numbers are hke the Bloch momenta in a 
periodic potential and label representations of the homotopy group of gauge 
transformations. On the fundamental domain the non-contractible loops arise 
through identifications of boundary points. Although slightly more hidden, 
the fundamental domain will therefore contain all the information relevant for 
the topological quantum numbers. Sufficient knowledge of the boundary iden- 
tifications will allow for an efficient and natural projection on the various su- 
perselection sectors. Typically we integrate out the high-energy modes, being 
left with the low-energy modes whose dynamics is determined by an effective 
hamiltonian defined on the fundamental domain (restricted to these low-energy 
modes) . In this it is assumed that the contributions of the high-energy modes 
can be dealt with perturbatively, generating the running coupling and the ef- 
fective interactions of the low-energy modes. We will in detail discuss the 
results for finite volumes with a torus and sphere geometry. 

3 Gauge Fields on the three- Torus 

Probably the most simple example to illustrate the relevance of the funda- 
mental domain is provided by gauge fields on the torus in the abelian zero- 
momentum sector. Let us first take G=SU(2) and Aj = i^Ta {L is the size of 
the torus, Tj are the Pauli matrices). These modes are dynamically motivated 
as they form the set of gauge fields on which the classical potential vanishes. It 
is called the vacuum valley (sometimes also referred to as toron valley) and one 
can attempt to perform a Born-Oppenheimer-like approximation for deriving 
an effective hamiltonian in terms of these "slow" degrees of freedom. To find 
the Gribov horizon, one easily verifies that the part of the spectrum for FP{A) 
that depends on C, is given by A^''(C) — 27rn • (27rn ± C)/i^, with n 7^ an 
integer vector. The lowest eigenvalue vanishes if Cj = ±2tt. The Gribov region 
is therefore a cube with sides of length 47r, centered at the origin, specified by 
|Cj | < 27r for aU j, see Fig. | 

The gauge transformation — exp{TrixjT3/ L) maps Cj to Cj + 27r, 
leaving the other components Ci^j untouched. As /iq) is anti-periodic it j.s 
homotopically non-trivial (they are 't Hooft's twisted gauge transformationsEI) 
We thus see explicitly that gauge copies occur inside fl, but furthermore the 
naive vacuum A = has (many) gauge copies under these shifts of 27r that 
lie on the Gribov horizon. It can actually be shown in the Coulomb gauge 
that for any three-manifold, any Gribov copy by a homotopically non-trivial 
gauge transformation of |4- = will have a non-trivial zero eigenvalue for the 
Faddeev-Popov operatorBJ Taking the symmetry under homotopically non- 
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trivial gauge transformations properly into aceount is crucial for describing the 
non-perturbative dynamics and one sees that the singularity of the hamiltonian 
at Gribov copies oi A = 0, where the wave functionals are in a sense maximal, 
could form a severe obstacle in obtaining reliable results. 

To find the boundary of the fundamental domain we note that the gauge 
copies C = (tt, C2,C3) and C = (— tt, C2,C3) have equal norm. The bound- 
ary of the fundamental domain, restricted to the vacuum valley formed by 
the abelian zero-momentum gauge fields therefore occurs where \C'j \ = tt, well 
inside the Gribov region, see Fig. ^. The boundary identifications are by the ho- 
motopically non-trivial gauge transformations . The fundamental domain, 
described by \Cj\ < tt with all boundary points regular, has the topology of 
a torus. To be more precise, as the remnant of the constant gauge transfor- 
mations (the Weyl group) changes C to — C, the fundamental domain A/G 
restricted to the abelian constant modes is the orbifold /Z2. As we will see, 
the orbifold singularities are associated with problems in using the harmonic 
approximation in perturbation theory for modes orthogonal to those of vacuum 
valley. 




Figure 3: A two dimensional slice of the vacuum valley along the (C\,C2) plane. The fat 
square give the Gribov horizon, the grey square is the fundamental domain. The dots at the 
Gribov horizon are Gribov copies of the origin. 



Formulating the hamiltonian on A, with the boundary identifications im- 
plied by the gauge transformations , avoids the singularities at the Gribov 
copies of A = 0. "Bloch momenta" associated with the 27r shift, implemented 
by the non-trivial homotopy of label 't Hooft's electric fiux quantum 
numbersO '^{Cj = — tt) = exp(7riej)^'(Cj — tt). Note that the phase factor is 
not arbitrary, but ±1. This is because /i? is homotopically trivial. In other 
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words, the homotopy group of these anti-periodic gauge transformations is Zf. 
Considering a shce of A can obscure some of the topological features. A loop 
that winds around the slice twice is contractible in A as soon as it is allowed 
to leave the slice. Indeed including the lowest modes transversa to this slice 
will make the Z2 nature of the relevant homotopy group evident E3 It should be 
mentioned that for the torus in the presence of fields in the fundamental rep- 
resentation (quarks), only periodic gauge transformations are allowedEj This 



will be discussed in Sec. 3.2 



In weak coupling Liischer showed unambiguously that the wave functionals 
are localiziejd around A = 0, that they are normalizable and that the spectrum 
is discretep'B In this limit the spectrum is insensitive to the boundary identifi- 
cations (giving rise to a degeneracy in the topological quantum numbers). This 
is manifested by a vanishing electric flux energy, defined by the difference in 
energy of a state with |e| = 1 and the vacuum state with e = 0. Although there 
is no classical potential barrier to achieve this suppression, it comes about by 
a quantum induced barrier, in strength down by two powers of the coupling 
constant. 

y(C) = 2L- + C| - |2.p|) . ^ 5: (8) 

The second identity is obtained by Poisson resummation@ and the periodicity 
reflects the gauge syiainetry discussed above. This quantum induced barrier 
leads to a suppressiono with|-a.iactor exp(— 5/5) instead of the usual factor of 
exp(— 87r^/g^) fniiJnstantonsEilS The associated tunneling configuration was 
called a pinchoncB At stronger coupling the wave functiopal spreads out over 
the vacuum valley and drastically changes the spectrumE^ At this point the 
energy of electric fiux suddenly switches on. 

Integrating out thejQon-zero momentum modes, for^which Bloch degener- 
ate perturbation theorycil provides a rigorous frameworkQ one finds an effective 
hamiltonian. Near A — 0, due to the quartic nature of the potential energy 
V{A) = — ^ / dPxtr {Fj^j) for the zero- momentum modes (the derivatives van- 
ish, such that the field strength is quadratic in the field), there is no separation 
in time scales between the abelian and non-abelian modes. Away from A = Q 
one could further reduce the dynamics to one along the vacuum valley, but near 
the origin the adiabatic approximation breaks down. This gives rise to a conic 
singularity, V{C) = 2|C|/i + O (C^), due to fluctuations of the non-abelian 
zero-momentum modes (contributing the p = term in Eq. (^)). 

The effective hamiltonian is expressed in terms of the coordinates c°, where 
i = {1,2,3} is the spatial index (eg = 0) and a = {1,2,3} the SU(2)-color 
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index. It will be convenient to introduce /"j, = —eabd.c'jcf and rj = \J^a '^j'^j- 
The latter are gauge-invariant "radial" coordinates, that will play a crucial role 
in specifying the boundary conditions. The zero-momentum gauge field and 
field strength read Aj{x.) = ic°j^ and Fjk = ifjk'it^- Fo^' dimensional reasons 
the effective hamiltonian is proportional to 1/L. It will furthermore depend 
on L through the renormalized coupling constant {g{L)) at the scale fi = 1/L. 
To one loop order (for small L) g^{L) = 127r^/[— 11 In(AMS'i)]- One expresses 
the masses and the size of the finite volume in dimensionless quantities, like 
mass-ratios and the parameter z = mL. In this way, the explicit dependence 
of (7 on L is irrelevant. This is also the preferred way of comparing results 
obtained within different regularization schemes (e.g. dimensional and lattice 
regularization). The effective hamiltonian is now given by 
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+72 ^"^t+l^^ ^1^] + 74 X! + 75 X! ''^^^l + 76 n "^'^ 

i i>j i i^j i 

+«3 E ^'(/j^)' +«4E^'(/5)' +«5det^c (9) 

ijk.a ij , a 

i i^j i>j i 

We have organized the terms according to the importance of their con- 
tributions. The first^two terms give (when ignoring ai^2) the lowest order 
effective hamiltonianP whose energy eigenvalues are O (g^^^) , as can be seen 
by rescaling c with g^^^. Thus, in a perturbative expansion c — 0[g^^^). 
Liischer's computationB to O (<?^^^), determined by ai,a2, 71,72 and 73, was 
extended to O (5^) , to ensure that the terms parametrized by % included the 
vacuum-valley effective potential (i.e. the part that does nO|Lvanish on the 
set of abelian configurations) to sufficient numerical accuracyts The order 
terms were just included to check stability of the results. The vacuum-valley 
effective potential could actually be computed at two loops to all orders in the 
field, up to an irrelevant constant 

T^2-ioop(C) = ^(AF(C) + AF(0)f, Ay(C)=E^^' (10) 

i 

but more cumbersome to compute were terms of the form g'^c^d^, also of two 
loop orderEa Their effect on the spectrum can, however, be neglected and 
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we have not listed these terms. The coefScients appearing in H^s have the 
following numerical values 



71 = 


-3.0104661 


10" 


-1 


- 3.0104661 


10" 




ai = 


+2.1810429 


10-2, 


72 = 


-1.4488847 


10" 


-3 


- 9.9096768 


10" 




a2 — 


+7.5714590 


10-3, 


73 = 


+1.2790086 


10" 


-2 


+ 3.6765224 


10" 


-'(5/2-)^ 


as = 


+1.1130266 


10-4, 


74 = 


+4.9676959 


10" 


-5 


+ 5.2925358 


10" 


-'(5/2^)^ 
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-2.1475176 


10-4, 


75 = 


-5.5172502 


10" 


-5 


+ 1.8496841 


10" 


-'(5/2^)^ 


as = 


-1.2775652 


10-3, 


76 = 


-1.2423581 


10" 


-3 


- 5.7110724 


10" 


-^(5/2-)^ 








77 = 


-9.8738947 


10" 


-7 


- 5.1311245 


10" 


-'(5/2^)^ 








78 = 


+9.1911536 


10" 


-6 


+ 9.1452409 


10" 


-'(5/27r)^ 








79 = 


-2.7911565 


10" 


-5 


- 2.5203366 


10" 


-'(5/27r)^ 
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+1.8208802 


10" 


-5 


+ 6.0939067 


10" 


''(5/27r)^ 
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A challenge was to rigorously determine the semiclassical expansion for the 
energy of electric flux due to the tunneling through a quaidum induced poten- 
tial barrier. The so-called Path DecompositioiLELxpansionpj was an important 
tool that could be tailored to this situationl£3'C3 One required matching the 
perturbative tail of the wave function derived from Liischer's effective hamil- 
tonian to the semiclassical contribution in the classically forbidden region. The 
resulting asymptotic expansion for the enepgiz-pf one-unit (|e| = 1) of electric 
flux, AE = Eoie) - EoiO) was found to beBB 



AE = 2L-iAB2g5/3(L)exp 



S/g{L) + s,T/g'^'{L) {l + /(g(L))}. (12) 



Here S — 12.4637 is the tunneling action and T = 3.9186 is related to the 
tunneling time, and comes from the classical turning points of the quantum 
induced potential V{C), which also determines A = 0.692.7, due to its trans- 
verse fluctuations along the tunneling path. It was showita that f{g) = (g). 
The constant Si = 4.116719735 determines to lowest order the mass of the 
scalar gluebaJL toq = L^^eig'^/^{L)^ which was already computed by Liischer 
and MiinsterEiJ Finally, the quantity B ~ 0.206 is taken from the asymptotic 
expansion of the lowest order perturbative wave function (in the direction of 
the vacuum valley). To compute ei and B it suffices to consider H^s in lowest 
non-trivial order, by putting all constants a; and 7; equal to zero. 

Of more practical importance is to go to the domain where the typical 
energies are above those of the quantum induced potential barrier, when the 
wave functional has spread out over the full fundamental domain, see Fig. ^ 
In this case one will become sensitive to the boundary identifications on the 



11 



fundamental domain. The influence of the boundary conditions on the low- 
lying glueball states is felt as soon as z = moL > 0.9. Wc summarize below 
the ingredients that enter the calculation. 

The choice of boundary conditions, associated with each of the irreducible 
representations of the cubic group 0(3, 'Z) and the electric flux quantum 
numbersp is best described by observing that the cubic group is the semidirect 
product of the group of coordinate permutations 5*3 and the group of coor- 
dinate reflections Z|. We denote the parity under the coordinate reflection 
c° —^ — c° (Va) by Pi = ±1. The electric flux quantum number for the same 
direction will be denoted by = ±1. This is related to the more usual ad- 
ditive (mod 2) quantum number ej by qj — exp(j7rej). Note that for SU(2) 
electric flux is invariant under coordinate reflections. If not all of the elec- 
tric fluxes are identical, the cubic group is broken to 52 x Zf, where 52 (= Z2) 
corresponds to interchanging the two directions with identical electric flux (un- 
equal to the other electric flux). If all the electric fluxes are equal, the wave 
functions are irreducible representations of the cubic group. These are the 
four singlets A^^y which are completely ( ant i-) symmetric with respect to 53 
and have pi — p2 — P3 ~ H. Then there are two doublets , also with 
Pi — P2 — P3 = ±1 and finally one has four triplets T^^y Each of these 
triplet sta±£s_can be decomposed into eigenstates of the coordinate reflections. 
Explicitlyli3'L3 for 7{(2) have one state that is (anti-)symmetric under in- 
terchanging the two- and three-directions, with pi = —p2 — —ps = ±1. The 
other two states are obtained through cyclic permutation of the coordinates. 
Thus, any eigenfunction of the effective hamiltonian with specific electric flux 
quantum numbers qi can be chosen to be an eigenstate of the parity operators 
Pi. The boundary conditions of these eigenfunctions 5'q,p(c) are given by 

*q,p(c)L,=, = , if piq, = -1 

^(n^'q,p(c))L^_ = , if p,q^ = +l (13) 

and one easily shows that with these boundary conditions the hamiltonian is 
hermitian with respect to the inner product < >= <^d9c^'*(c)^''(c). 



In Sec. 3.1, when discussing the case of SU(3), we show in more detail how one 
arrives at these boundary conditions. Fon_negative parity states (Jli-P* = ~1) 
this description is, however, not accurateEj as parity restricted to the vacuum 
valley is equivalent to a Weyl refiection (a remnant of the invariance under 
constant gauge transformations). One can use as a basis for ^'p_q(c) 

3 

< c|l, n; e >= ^ W{\, m)Y[r-'xti{n)Yu,r.M, ^.), (14) 
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with Yi^m spherical harmonics, and W{1, m) the Wigner coefficients for adding 
three angular momenta X)i = ^ (ensuring gauge invariance). The angular 
quantum numbers k are restricted to be even or odd, for resp. pi — 1 or 
Pi — —1, and the radial wave function X^'^)(r) either vanishes, or has vanishing 
derivative a,tr — n (for resp. l + e even or odd, see Eq. (p^)). One computes the 
matrix elements of the effective hamiltonian (Eq. (^)) and solves for the low- 
lying spectrum by Rayleigh-Ritz (providing also lower bounds from the second 
moment of the hamiltonian). The first two lines in Eq. (|^) are sufficient to 
obtain the mass-ratios to an accuracy of better than 5%. 



L.V(A) 




Figure 4: Sketch of the potential for the torus. Shown are two vacuum valleys, related to 
each other by a gauge transformation hi, with winding number u{hi) = 1. The induced 
one-loop effective potential, of height 3.210/L, has degenerate minima related to each other 
by the anti-periodic gauge transformations The classical barrier, separating the two 

valleys, has a height 72.605/Lg^. 

At larger volumes extra degrees of freedom will behave non-perturbatively. 
We know from the existence of gauge transformations with non-trivial winding 
number 

y{h) = ^ j^^ir{{h-Hhf), (15) 

that there exist gauge equivalent vacuum valleys, that can be reached only by 
crossing a saddle point (called the finite volume sphaleron), see Fig. ^. Typ- 
ically this saddle point lies on the tunneling path (the instanton), with the 
euclidean time of the instanton solution playing the role of the path parame- 
ter. We expect, as will be shown for the three-sphere, that the boundary of 
the fundamental domain along this path in field space across the barrier occurs 
at the saddle point in between the two minima. The degrees of freedom along 
this tunneling path go outside of the space of zero-momentum gauge fields 
and if the energy of a state flows over the barrier, its wave functional will no 
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longer be exponentially suppressed below the barrier and will in particular be 
non-negligible at the boundary of the fundamental domain. Boundary identifi- 
cations in this direction of field space now become dynamically important too. 
The relevant "Bloch momentum" is in this case the parameter. Wave func- 
tionals pick up a phase factor e'^ under a gauge transformation with winding 
number one. 




0.5 1 1.5 2 2.5 3^ 0.5 1 1.5 2 2.5 3.. 

->9 ->£r 

Figure 5: The top figures show LE{g) for the relevant (positive parity) representations. The 
dotted hnes are without the two loop correction included. The dashed curve denotes the 
barrier height LE^pfi{g). The bottom figures show z(g) = Lm(g) = L(E{g) — Eo{g)), and 
the dashed line is at z = 2it. 

To estimate for which volumes the extra degrees of freedom start to con- 
tribute non-perturbatively, the minimal barrier height that separates two vac- 
uum valleys was found to be Egphig) — 72.605/Iig^, using the lattice approxi- 
mation and carefully taking the continuum limitEj As long as the states under 
consideration have energies below this value, the transitions over this barrier 
can be neglected (or treated semi-classically if there is no perturbative con- 
tribution) and the zero-momentum effective hamiltonian provides an accurate 
description. One can now find for which volume the energy starts to be of the 
order of this barrier height. For this purpose we have collect the energy levels 
as a function of g in Fig. |^. In the top two figures we have plotted LE{g) for 
the various representations of interest (0 stands for the vacuum state, which 
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is in the representation of the cubic group). The dotted curves are with- 
out the two loop corrections included. In the same figures the dashed curve 
shows the height of the instanton barrier. In the lower two figures we show 
z = Lrn = L{E — Eq). The dashed line is at z = 27r. We see that the sensitivity 
to the two loop corrections almost entirely drops out for the energy differences, 
at the scale of this figure. We now read off that instantons only become im- 
portant for zq = 5 to 6, i.e. for L roughly 5 to 6 times the correlation length 
set by the scalar glueball mass. Certain states may actually be less sensitive to 
instantons, in case the representation is such that the wave functional in the 
direction of the barrier is suppressed. This may for example be the case for the 

state. It should also be noted that z = 2tt is associated with the energy 
scale 27r/L of the non-zero momentum field modes that have been integrated 
out. Therefore, z should also not be much bigger than 27r, although also here 
the particular representation may matter. 

On the three-torus we have therefore achieved a self-contained picture of 
the low-lying glueball spectrum in intermediate volumes from first principles 
with no free parameters, apart from the overall scale. For very small volumes 
the energy of electric flux vanish and there is an accidental rotational symme- 
try, only split by the O (g^^^) terms in Eq. (^). The 2+ tensor glueball splits 
in a nearly degenerate doublet and a triplet T^. Both are lighter than 
the scalar 0+, also denoted by as the scalar singlet representation of the 
cubic group. Going to larger volumes, L > 0.1 fermi, the energy of electric 
flux per unit length, which in an infinite volume would be the string tension 
K, is surprisingly constant in intermediate volumes, whereas the splitting of 
the tensor states becomes quite large. In intermediate volumes the doublet 

and triplet T2 have respectively masses of roughly 0.9 and 1.7 times the 
scalar mass. This doublet E~^, lighter than the scalar, was first observed in the 
lattice studies of Berg and Billoireo and caused some stir at the time. 

The lattice data for-the triplelH were obtained only after our continuum 
results first appearecCJ'Ea and did not confirm the predictions for this state. 
This was resolved by Vohwinke£3 by observing that the state we had initially 
identified as (taking pi — P2 = P3 = ^ instead of pi ~ —p2 = —pz = 1) 
actually carried two units of electric flux (cmp. Eq. (|l3|)), making it even 
more of a surprise that it was found to be even lighter than the doublet E^ 
in intermediate volumes. In the infinite volume limit it is pushed out of the 
low-lying spectrum. Around the same time this state (named Tn) was as such 
measured on the latticepi confirming the proper interpretation of these states. 

Electric fiux energies (for the trivial representation) are labelled by e = ei 
for e = (1, 0, 0), 62 for (1, 1, 0) and 63 for (1, 1, 1). For e„ we speak of n units of 
electric flux, sometimes called "torelon" energies. As was argued by 't Hooftp 
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if a confining string would have formed, = KL (or = LE^ = KL^), it 
would be energetically favorable to run along the direction of e, giving i?„ = 
Eicn)/ E{ei) = instead of splitting in separate strings, each winding in 
the direction for which Cj — 1, which would gpe i?„ = n. In intermediate 
volumea.it is the latter behavior that we foundpl confirmed by Monte Carlo 
resultsE3 Onedias to go to quite large volumes to start to observe the expected 
^/n behaviorLj The same is found for SU(3), in a study where ajiifferent lattice 
Monte Carlo method was used to get the electric flux energies^ 




1 2 3 4 5 0.5 1 



Figure 6: Left: The ratio ^ K{L)/mE as a function of ze = mE(L)L. Full (contin- 
uum), dashed-dotted (6^ lattice) and dashed (4'^ lattice) curves give the hamiltonian results. 
Rayleigh-Ritz diagonalization errors are smaller than the thickness of the lines. The sensi- 
tivity to the two loop correction can be read off from the dotted curve (4^ lattice), where 
this correction was not included. Note the bloYgrW scale. Right: electric flux ratios R„ 
as a function of Ze = Ee{L)L. Monte Carlo dataOCj obtained with the Wilson action on a 
Af3 X Aft lattice {Nt > 6N), for N = i (dots), 6 (squares), 8 (triangles) and 10 (diamond). 

The lattice Monte Carlo data of Berg and Billoir(S@ are compared with 
the hamiltonian results in Fig. ^. On the left we show the ratio \J K(L) j ze 
as a function oi ze- Their methods had considerable difficulty in dealing 
with the scalar glueball state, since it has the same quantum numbers as the 
vacuum. To push the results to low values of z^;, they used lattices N'^ x Nt 
with Nt = 256 for iV = 4. Please note that z — tueL around 0.95 is nearly 
a constant (even not single valued) function of g{L), see Fig. |[ This over- 
emphasizes the steep behavior where the wave function starts to spread over 
the whole vacuum valley, leading to non-zero energies of electric flux. We took 
the Monte Carlo data from Table V of their paperpi but removed those entries 
with Nt/N < 6. Qa the right in Fig. || we compare with the lattice Monte 
Carlo data of BergcZl for the electric flux ratios, but used where available, the 
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more accuxate resultalj for = LEe- The data at/3 = 4/5Q = 2.7 from both 
papersEIrE3 were left out because they were not consistent with each other. 

Considerable progress was achieved by the so-called variational method, 
which allowed much more accurate results, not only for the scalar glueball, but 
also for arbitrary other representationsE^I Jn Fig. |^ we present a comparison 
with the Monte Carlo results of Michaelpll obtained for^ lattice of spatial 
size 4'^, confirmed in a study of improved lattice actions 13 The hamiltonian 
results below z = 0.95 are due to Liischer and Miinster^ which is where the 
spectrum is insensitive to any identifications at the boundary of A. The electric 
flux energy ratios i?2 and i?3 are shown on the right, cmp. Fig. ^. 
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Figure 7: Left: mass ratios m/mo as a function oi z = mo{L)L, with mo the scalar glueball 
mass. Full (continuum) and dashed (4"^ lattice) curves give the hamiltonian results. Shown 
are string tension yj' K{L), and tensor masses and the exotic states Tn with two, and 
Till = 72(111) with thnen units of electric flux. Right: electric flux ratios ij„ as a function 
of 2e. Monte Carlo dataEj obtained with Wilson action on 4^ X Nt lattice {Nt = 32, /3 = 2.4 
and Nt = 99, /3 = 3.0). 



We note that the solid curves that represent the continuum results, which 
were reproduced by Berg and Vohwinkelpio deviate significantly from the 
lattice data. Initially the lattice data were not accurate enough to show this 
deviation. Even though one should not expect a 4^^ lattice to be an accurate 
approximation for the continuum, it was cause for some doubt that the ap- 
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proximations made in the eontinuum studies were not under control as wellEI 
To settle this issue we redid the complete derivation of the effective harniJixtj 
nian starting from the lattice theory, without taking the continuum lim,itL5'cB 
The hamiltonian is basically of the same form as in Eq. (^) , except that the 
coefficients in Eq. (pi] ) depend on the lattice spacing (some extra corrections 
appear, e.g. to correct for the discrete time evolution). One can follow the 
renormalization group flow of the hamiltonian to its continuum fixed point in 
this formalism in all detail, see Fig. ||. Using the same analysis as in the con- 
tinuum leads for a finite lattice to the dashed curves in Figs. 0,^. The lattice 
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Figure 8: Flow of the lattice hamiltonian coefficients to their continuum (A^ — + oo) values. 
Shown are ai{N)/ai (dashed lines) and '~ii(N)/'~fi (full lines) as a function of the square of 
the lattice spacing (afV)^ = 



data now agree perfectly, up to a volume of about 0.75 fermi, the regime for 
which we have shown that the effective hamiltonian in the zero-momentum 
modes should provide a good approximation. The deviations for i?2 and in 
particular for are in accordance with the fact that these are of relatively 
high energy, and therefore expected to be sensitive to other non-perturbative 
effects that invalidate integrating out all the non-zero momentum modes, see 
Fig. |. 



3.1 General Gauge Group 

The question of extending the previous results l(mSU(3) is a natural one 
in the ligii±_,of QCD. The perturbative expansion|jL3 vacuum-valley effective 
potentiaiP'Ej and the semi-classical evaluation of the energy of electric fluxo 
(due to the tunneling through the quantum induced vacuum-valley effective 
potential) are more or less straightforward. The non-trivial problem of for- 
mulating the appropriate boundary conditions on the boundary of the funda- 
mental domain for SU(3) was solved by Vohwinkel and qualitative agreement 
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with the lattice data was foundH ci-With the results for SU(3) in hand gen- 
eralization to SU(Af) was achievedc3 The formalism sketched in Sec. |^ was 
not yet developed in those early days. What is now called the fundamental 
domain was then called the unit cell. Although the argumentation was more 
cumbersome, the underlying principles were the same. 

A simple approximation for Liischer's effective hamiltonian in terms of the 
zero-momentum gauge fields Aj{x.) — ic'^Ta/L {Ta the hermitian generators of 
SU(iV) and fij = L'^Fij{x) the dimensionless field strength) is given by 

^e«(c) ^ E - + -) E (/J) + V.icl (16) 

with 

Viic) = I E (v/Trad(2^P + c''adr,)2- Tr ad(2^p)2) , (17) 

where g is the renormalized coupling constant at the scale ji = 1/L. In a 
perturbative expansionPo this gives the correct result up to O {g^^"^)- Along 
the vacuum valley, parametrized by Aj{x.) — i^^Si^ C^Th/L, where the first 
A'^ — 1 generators are assumed to commute (forming a basis of the Cartan 
subalgebra He), the effective potential Vi{c) is exact to one loop order. A 
more explicit result can be found in Sec. |3.2| , Eq. (p5|). As for SU(2), one does 
not include the p = term, since only the non-zero momentum modes are to 
be integrated out. 

The effective potential Vi{c) only depends on the Casimir invariants 

^2tr[{ctTaf], ^ Atr [{c^Taf] , (18) 

where the sum over all color indices is implicit. There are as many inde- 
pendent Casimir invariants as the rank of the gauge group. For SU(2) only 
ri will be non-trivial. These coordinates are uniquely related to those ob- 
tained by restricting the zero- momentum gauge field to be abelian. For SU(2), 
c^Ta = C,T3/2 yields rf = Cf, whereas for SU(3), in terms of the Gell-Mann 
matrices, = (C^As + Cf\3)/2 gives s, = ClVSiiCf)^ - (01)^3) and 

rj = (Cj^)^ -I- (Cf)'^. In this manner the effective potential on the set of abelian 
zero-momentum modes can indeed be minimally extended in a unique way to 
all constant gauge fields. 

By adding the zero-momentum (p — 0) contribution to the one-loop effec- 
tive potential, 

Kff(c) = V^i(c) + 2L-VTrad[(c'^ra)2] , (19) 
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restricted to the vacuum valley gives the appropriate effective potential when 
integrating out all the field modes orthogonal to the vacuum valley. Its sym- 
metries are a consequence of gauge invariance, which can be divided in two 
classes. The constant gauge transformations, that leave Hq invariant. This is 
represented by the Weyl group W acting on C** (for SU(2) C — C). The 
other class of gauge transformations that leaves the set C** invariant are of 
the form /i0(x) = exp(27rix • @/L), where 8i G Hq such that h does not 
affect the periodic boundary conditions on the gauge fields. These lead to shift 
symmetries on C**. The associated lattice of Qj corresponds to the dual root 
lattice A,., as follows from the condition that exp(27ri0j) = 1. The vacmmi 
valley, i.e. the moduh space of flat connections {Fij = 0), corresponds to the 
orbifold [Hg / 2tt Ar]^ /W. Recently it has become clear that for orthogonal 
and exceptional Lie-groups there are other, disconnected, components in the 



moduli space of flat connections, which we will briefly discuss in Sec. 3.5. 

A crucial role is played by the twisted gauge transformations, only periodic 
up to an element of the center of the gauge groupB These can be realized 
with h&, but now Qj belongs to the dual weight lattice A^, which follows from 
the requirement that exp(27rj0j) is an element of the center Zg of the gauge 
group. Indeed, A^j/Ar = Zg- These twisted gauge transformations generate 
a shift symmetry associated with A^. Their homotopy type is specified by 
G (A^/Ar)^ 9i Zq. A particular representative for SU(A'^) is given by 
/ik(x) = exp(27rik • xQq/L), where Oq G Hg is a generator for the center, 
exp(27rz6o) = ex.p{2TTi / N) . For SU(2) we can choose 9o = T3/2 and for 
SU(3) there are two independent choices (that span the dual weight lattice A^,) 
Go ^ Ag/VSande^, = ^(Aa-As/VS). The non-trivial homotopy is labelled by 
k G Zp^ and is thus in general non-trivially reppsented on the wave functional. 
Any gauge transformation can be decomposedE2l in h = hkhi{h)hQ, where hi is 
a particular strictly periodic gauge transformation with unit winding number, 
I' (hi) — 1. What is left is a homotopically trivial gauge function ho, under 
which the wave functional is invariant. We therefore have 

•^i[h]A) ^ *([/ik/ii]^) = exp(27rie • k/N) expiiOiy)^ (A) , (20) 

where 6 is the usual vacuum parameter associated with instantona and e 
(defined modulo N) is the gauge invariant definition of electric fiuxlj As for 
SU(2) the electric flux quantum number can be implemented within the zero- 
momentum effective hamiltonian by imposing suitable boundary conditions on 
the boundary of the fundamental domain. For this we study the action of 
/ik on C^ corresponding to a shift. For SU(2) C ^ C + 27rk and for SU(3) 
(Ci , C2) ^ (Ci -h 27r(2n - 1)/ V3, + 27rl) , with k = n + 1 (mod 3) specifying 
the homotopy type of the gauge transformation. 
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Fig. ^ gives for SU(3) the equipotential lines of the effective potential VcS, 
in a plane specified by putting C2 = C| = 0, the fundamental domain A 
(bounded by the dashed hexagon) and the Gribov region f2 (bounded by the 
fat hexagon). The triangular lattice structure is due to the invariance under 
twisted gauge transformations and reflects the dual weight lattice of SU(3). 
Indeed, the gauge transformations exp(2TriQoXj / L) and exp{2TriQ'QXj / L) all 
map points on the boundary of the fundamental domain to the same boundary. 
They also map = to the Gribov horizon. The Weyl transformations are 
generated by the reflections in the three principle axes of the dual weight 
lattice. 




Figure 9: Two dimensional cross- section of V^ft, for C^T,, = A8/2 + A3/2 and C| 3 = 0. 
The circles indicate the equipotential lines (in increments of Vcff{AiT5)/7). The light shaded 
area represents the fundamental region, with the dashed hexagon as its boundary, the fat 
hexagon indicates the Gribov horizon, with gauge copies of yl = indicated by the dots. 
The darker shaded area is the so-called fundamental Weyl chamber. 

The near spherical behavior of the equipotential lines within each triangle 
is due to a remarkably accurate approximation (cmp. Fig. that holds for 
any SU(A^) when restricting Cf to one dimension, say Ci = (normalizing 
as usual tr(raTf,) = ^Sab) 

Kff(C) « N {{AnSf - (C - AnSf) /2^i, (21) 

for C in the Weyl chamber restricted to the fundamental domain, centered 
around 47r(5, where 6 is the highest weight. For SU(2) 5 = 1/4 and for SU(3) 



21 



S = {6^,S^) = (\/3, l)/6. It is also instructive from Fig. ^ to identify the 
orbifold singularities, along the axes of the dual weight lattice (edges of the 
Weyl chambers), which are fixed under the Weyl reflections. This is where 
the U{1)'^ symmetry in a generic point of the vacuum valley is restored to 
SU(2), whereas at the corners (all gauge copies of A = 0) the full SU(3) 
symmetry is restored. At these locations also some of the directions in which 
the classical potential is quadratic turn quartic. It is this that gives rise to the 
conic singularities in Vcs- The effective potential Vi in Eq. ([l7|), when restricted 
to the fundamental domain, is free from any of these singularities, since they 
are caused by fluctuations in the non-abclian constant modes, which are not 
integrated out in Eq . (p7|). The con ic singularities of VcS in A are therefore all 
described by 2L~^ Tr ^d[{C''Tb)'^] = 2NL-'^\SbC''\, see Eq. (||). 

To implement the symmetries on the wave functional we perform a change 
of coordinates cf (Cf,r2i), where fii stands for the collection of SU(A^)- 
angular coordinates parametrizing SU(7V)/U(1)^~-^ and (Cf) are restricted 
to a fundamental Weyl chamber, for SU(3) the triangular shaded region in 
Fig. ^, where the relation between (C^, C^) and (r, s) is one to one. For SU(2) 
these new coordinates are the spherical coordinates (r^, fi^), with fli = {6i, (pi) 
coordinates on = SU{2)/U{1). The square root of the Jacobian of this 
transformation, is given by J = J^j Ji, with J,j = Q for SU(2) and = 
Cf{{Cfy — 3(C/)^) for SU(3). The wave functions can be decomposed as 

^{c)^l[jr\cr)xAC^)ydn,). (22) 

i 

The "radial" part x will be antisymmetric with respect to Weyl reflections, so 
as to cancel the zero's of the Jacobian. For SU(2) the angular wave functions 
y{fl) are nothing but the spherical harmonics, see Eq. (p^, whereas in general 
they are irreducible representations of the gauge group. Helpful in making suit- 
able choices for x is that the vacuum valley kinetic term after the rescaling with 
the square root J of the Jacobian becomes —^g'^^i^^d/dCfY, compatible 
with the canonical flat metric on the orbifold. Furthermore, as is familiar from 
the radial reduction in three dimensions, this rescaling does not create an ad- 
ditional potential term, since J is in all cases harmonic, ^ J / {dCf Y — 0. 

We note that suitable combinations of a shift and Weyl reflection leave 
invariant the lines that constitute the polygon at the boundary of the funda- 
mental domain, see Figs. ^ and ^. This implies that alternatively the proper- 
ties of X can be described in terms of boundary conditions at the boundary 
of the fundamental domain. E.g., for SU(2) this leads to two possible choices 
of boundary conditions at Vj = tt, for each j, see Eq. (^). One can use 
x(C) = sin(nC/2), although spherical Bessel functions provide a more effi- 
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cient choice, keeping the hamiltonian more spaxseJp CarefuUy working out the 
consequences of the symmetries one can showE§CJ that a complete basis for x 
in the case of SU(3) is given by 



xiC^.C^) = sin(mCV2)exp(mCl^/3/2) (23) 
+ sin (m(\/3Ci - C^)/^) exp [in\/?,{ViC^ + C^)/2^ 
- sin (m(V3Ci + C^)/A^ exp [in^/2,{^/iC^ - C^)/2^ , 

for each of the three coordinate directions. The quantum numbers n and m will 
be restricted by the electric flux and irreducible representations of the Weyl 
and cubic groups, which is the part that requires most of the care. Finally 
if one restricts Yl - to transform as a singlet under SU(iV) one obtains a 
complete and gauge invariant basis for the effective hamiltonian, that through 
the boundary conditions carries the information of electric flux. Needless, to 
say that for SU(3) the computation of all the relevant matrix elementscS is 
rather more cumbersome than for SU(2). However, once the matrix of the 
hamiltonian for a particular basis is computed one performs a simple Rayleigh- 
Ritz analysis to determine the spectrum. The region where the wave functional 
spreads out over the whole vacuum valley, with the energy of electric flux no 
longer suppressed by the quantum induced barrier of Voff, is well described by 
this Rayleigh-Ritz analysis, based on Eq. (p3|). Of course, also for SU(3) the 
results are valid as long as the classical barrier is sufhciently high as compared 
to the energy-levels studied, cmp. Fig. ^. Like for SU(2) the approximations 
break down for L large than 5 to 6 times the correlation length of the scalar 
glueball. That beyond this volume, L 0.75 fermi, instanton effects sei-in can 
be seen from the rather sudden onset of the topological susceptibilityEl 

3.2 Including Massless Quarks 

Quarks fields are not invariant under the center of the gauge groups. This 
means that on the space of zero-momentum abelian gauge fields, which still 
form the classical vacuum valley, the twisted gauge transformations no longer 
represent a symmetry. This complicates matters since without the equivalence 
under twisted gauge transformations the fundamental domain extends to the 
Gribov horizon. Some of these complications may be avoided by introducing 
on the original fundamental domain additional quark fields, obtained by ap- 
plying a twisted gauge transformation. The boundary identification on the 
boundary of the fundamental domain are include an operation that permutes 
these fermion field components. This has not been worked out so far. 
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Nevertheless, interesting statements can be made about the vacuum struc- 
ture in small enough volumes, for which the wave functional is sufficiently lo- 
calized around the vacuum configuration. One simply adds in one loop order 
the quantum effects of the quark field fluctuations. The resulting effective po- 
tential will no longer respect the center symmetry, but it still properly reflects 
invariance under constant and periodic gauge transformations. The quark 
fields can satisfy either periodic or anti-periodic spatial boundary conditions. 
Actually, for SU(2) (with —1 a non-trivial clement of Z2) these are equivalent 
by a twisted gauge transformation with homotopy type k = (1, 1, 1). Under 
this gauge transformation the gauge field is shifted and shows it is not a priori 
clear that ^ = will represent the proper classical vacuum to expand around. 
As we will show, it will be the correct one only with anti-periodic boundary 
conditions for the quark fields, both for SU(2) and SU(3). In that case, due to 
the anti-periodicity, there will be no zero-energy modes for the quark fields, and 
chiral symmetry is unbroken in the finite volume. For SU(3) no gauge equiva- 
lence of periodic to anti-periodic boundary conditions holds, and the vacuum 
structure with periodic quark fields actually leads to spontaneous breaking of 
some discrete symmetries. Yet, no zero-energy quark modes appear, and chiral 
symmetry remains unbroken. It also means that in a small volume, with quark 
momenta of the order tt/L and glueball masses of order g'^/'^{L) / that glue- 
balls cannot decay in mesons. The quark degrees of freedom can be integrated 
out. 
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Figure 10: Vacuum-valley effentive potential at C2 = C3 = 2tt for SU(2) and two flavors of 
periodic massless quarks (Eq. ( psj ) with k = 0), normalized to at its minimum, as well as 
the effective potential for zero flavors (Eq. (js])). The dashed curve represents Eq. (|2l[). 

To be more specific let us first generalize the computation of the vacuum- 
valley effective potential to include the quark field fluctuations. The most 
efficient way to represent the result is to introduce the weight vectors /x''^ 
determined by the eigenvalues of the abelian generators, 

T, = diag(Mi'\A.i2\...,Mf'), neHG. (24) 

For SU(2) one finds /i^^^ = i(l, -1), whereas for SU(3) ^t^' = (1, 1, -2)/VT2 
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and /i^^) = 1,0) (the conventions used in this paper-axe Ti — iAg and 

T2 = ^^3)- With Uf flavors of massless quark fields we findtj (see Fig. |l^) 

KKc^) = - /^i'^]) - "/ E ^(cVi^^ + ^k), (25) 

i>j i 

with k = or k = (1, 1, 1), for rasp, periodic or anti-periodic boundary con- 
ditions on the quark fields. The function ^(C) is the SU(2) one-loop effective 
potential for n/ = 0, Eq. (||). The correct quantum vacuum is to be found at 
the minimum of this effective potential. 

Observe that the gauge symmetry should not be spontaneously broken, 
which implies that the Polyakov loop observables 



P, (x)^ltr (^^cxp(^ 



dsAj(x + se(j)) , (26) 



should be a constant center element at the vacuum configurations, or 

P, = ^ tr (exp(zC)Tfc)) = 1 ^ exp(z;.("'c)) = exp(27rzZ,/iV). (27) 

n 

This implies that /i["-'C'' = 2-k\/N (mod 27r), independent of n, and gives 
= — 7i/AfV^(27rl/iV-f 7rk). In the case of anti-periodic boundary conditions, 
k = (1,1,1), this is minimal only when 1 = (mod 2-k). This means the 
quantum vacuum in this case is the naive one, A — Q {Pj — 1). In the 
case of periodic boundary conditions, k = 0, the above candidate vacua have 
1 7^ 0, that is Pj correspond to non-trivial center elements. Both for SU(2) 
and SU(3) this means Ij = ±1. For SU(2) the vacuum with Pj = —1 is unique, 
as also follows from the gauge equivalence argument given above. The only 
difference is that one now needs to expand around C = 27r(l, 1, 1). For SU(3), 
however, there are now 8 possible choices Pj = exp(±27ri/3), related by the 
three coordinate reflections. As this is a symmetry of the full hamiltonian, 
each is indeed equivalent. But it does mean that in a small volume parity (P) 
and charge conjugation (C) are spontaneously broken, although CP is still a 
good symmetryE£l A consequence of the spontaneous breaking of parity is that 
the mass gap, the lowest excitation above the vacuum, is exponentially small 
in a small volume. All these intricate effects would make this an ideal testing 
ground for dynamical fermion algorithms in lattice gauge theory. 

The minima of the effective potential are obtained from A = by the 
twisted gauge transformation /ii(x). As there are no zero-energy quark fleld 
modes, also the effective hamiltonian can be expressed as in the bosonic case 
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in terms of the zero-momentum gauge fields, after taking tliis shift into ac- 
count. In the case of SU(3) with periodic boundary conditions, because of the 
spontaneous break down of parity and charge conjugation invariance, extra 
terms appear in Liischer's effective hamiltonian. To the order in which this 
hamihonian was worked out the new interaction tr ({Ai^ Aj}Ak) appears. In 
addition couphngs that were equal because they were related by the discrete 
symmetries, now spontaneously broken, will splitllj All these corrections are 
linear in the number (rif) of quark flavors. The flavor dependence of the ef- 
fective hamiltonian for the case of anti-periodic boundary is simpler, since no 
new couplings appear. For SU(2) one simply replaces and 7^ in Eq. by 
ai — 2rija^ and 7,; — 2nf^'^, with 
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(28) 

Independently Kripfganz and Michael calculated for SU(A^) to O {g^l^) the 
change in the coefficients of Liischer's effective hamiltonian, due to quarks 
with anti-periodic boundary conditions onlytSI They confirmed the values in 
Eq. (|2^) and also introduced for SU(2) a lagrangian formulation of the effective 
hamiltonian in terms of compact group variables, that incorporates in a simple 
way the propsiLboundary conditions in field spaceEj After full equivalence was 
establishedpia the Monte Carlo analysis of this effective lagrangian model 
continued tcj_suffer from a technical difficulty in efhciently implementing t 
kinetic termpi only fully understood by Vohwinkel a number of years later 
This has hampered using the lagrangian formulation as a reliable alternatively 
for the hamiltonian Rayleigh-Ritz analysis. 

There is another choice of boundary conditions that strongly reduces the 
center symmetry. This is called C periodic boundary conditions, where the 
field is periodic up to a charge conjugation. The boundary conditions can 
be used to avoid the system to be charge neutral, as is the case for periodic 
boundary conditions, both for magnetic and electric charge. This was the 
originaL motivation to introduce these boundary conditions for the abelian 
theoryp3 which were subsequently studied for non-abelian gauge theoriesEfl 
For SU(2), which is pseudo real, one retrieves the periodic case, but for SU(3) 
the center symmetry is completely broken and a number of the features we 
saw when including quark fields, appear here as well. The vacuum valley for 
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SU(3) is reduced frcpa six to three dimensional in terms of the real gauge field 
A = iiCAa/i, withka 

Kff(C) = F(iC) + ViC - nl) + y(iC - ^1), (29) 

where V{C) is the SU(2) effective potential, see Eq. (|), and 1 = (1, 1, 1). The 
minimum of this effective potential occurs at the four points C — (7r,7r,7r), 
(— TT, — TT, tt), (tt, — tt, — tt), and (— tt, tt, — tt), which correspond to orbifold sin- 
gularities with quartic modes (associated with the three real generators iXa/^, 
for A — 2, 5 and 7, forming an SU(2) subalgebra). The effective hamiltonian is 
again of the Liischer type, at O (g^^'^) identical to it, at higher order additional 
couplings appear because of the spontaneous breaking of parity, i.e. the cubic 
group is broken down to the permutation group. No attempts have been made 
to study the fundamental domain for this theory, and we may expect similar 
difficulties as in the presence of quark fields. 

3.3 The Renormalized Coupling 

We have assumed there is a renormalized coupling in terms of which per- 
turbation theory in the field modes that are integrated out is well-behaved. 
By expressing quantities in dimensionless combinations, lattice Monte Carlo 
results and continuum (or lattice) hamiltonian results can be compared with- 
out being sensitive to any problems in expressing the renormalized coupling in 
terms of the bare coupling. Determining the renormalized strong coupling con- 
stant non-perturbatively in a reliable way is, however, an important problem. 
The integration constant of the beta-function, the so-called Lambda parame- 
ter, ideally should be fixed in terms of the infrared quantities of the theory, 
like the mass gap and string tension, or other observables in the low-lying 
spectrum of the theory. The running of the coupling allows one to compute 
unambiguously the strong coupling constant at, say the Z-boson mass. The 
most acmrate such method is based on a finite volumeL_atudy, proposed by 
Liischer pi long before it was feasible to be implementedl£^ It makes use of a 
discrete version of the beta-function, the so-called step scaling function. The 
scale at which the renormalized coupling is defined is fixed by the volume. The 
volume is subsequently changed by an integer factor (usually, but not neces- 
sarily) of 2. So instead of an infinitesimal scale transformation it considers 
a finite one. The change in the coupling can of course be obtained by inter- 
grating the beta-function, but this function is not available non-perturbatively. 
Instead, one picks a suitable definition of the renormalized coupling constant 
at the scale L, set by a given physical volume, than doubles L and calculates 
the value of the coupling at this new scale 2L. This can all be performed on 
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the lattice (hence the integer scahng factor) using Monte Carlo calculations, 
at each step carefully extrapolating to the continuum limit. Also (euclidean) 
time is taken finite, T = L. Small volumes go together with high temperatures 
in such geometries. The naive strategy of taking L large, and defining an ob- 
servable set by a variable scale within the same finite volume, fails because the 
lattice spacing gives a limit to the shortest distance one can probe in a given 
volume, due to computer limitations. 

Many definitions of the renormalized coupling could in principle be used, 
but technical requirements have led to a particular one that is related to the 
effective action in the background of a|-eonstant chromo-electric field, based 
on the so-called Schrodinger functional^ (SF). In the spatial directions the 
boundary conditions are periodic, but in the (finite) time direction one pre- 
scribes an initiaLand a final configuration of gauge fields taking values in the 
vacuum valleyE^IIS 

SU(2) : Cj(t = 0) = 27r(5 + 27/, Cj{t = T) = 2TT-Cj{t = Q), 
SU(3): (C},C|)(t = 0) = 27r(5+(-V377,0), (30) 
(C}, C|)(t = T) = (4V^/3, 0) - {C],C^){t = 0), Vj, 

with 5 the highest weight, defined below Eq. (^T|). The classical equations of 
motion lead to a linear interpolation in time, giving rise to a constant chromo- 
electric background. The euclidean quantum effective action, r(77), describes 
the reaction of the system to this background. The renormalized coupling 
is defined by g^{L) — ^^^-^ / ^^^-^^ evaluated at 77 = 0, where ro(ry) is the 
bare effective action. This background has been chosen to stay well away 
from the orbifold singularities along the vacuum valley for the entire classical 
path that interpolates between the initial and final configuration, to simplify 
the perturbative analysis (in part for estimating the lattice artifacts). This 
particular choice of coupling fits extremely well to the perturbative running of 
the coupling constant up to the largest volume probedq with L up to ~0.35 
fermi. At large volumes the non-perturbative running is bound to deviate. 

An alternative definition for the non-pcrturbatively defined running cou- 
pling for SU(2) has been based on ratios of the expectation values of suit- 
able Polyakov loop operators, using twisted boundary conditionsp the so-called 
twisted Polyakov (TP) schemelEl We will see next that twisted boundary con- 
ditions remove the zero-momentum modes, making perturbation theory well 
behaved. Without these twisted boundary conditions, computing expectation 
values in the four dimensional euclidean finite volume is difficult to controlEj 
The TP coupling also agrees well with perturbation theory and after matching 
of the scales, with the SF couplingEa Only the largest volume result [L = 0.28 
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fermi) probed by the TP coupling lies slightly, but significantly, below the per- 
turbative result. The near perturbative behavior seems to support the fact that 
non-zero momentum modes do behave perturbatively in intermediate volumes. 



3.4 Twisted Boundary Conditions 

With twisted boundary conditions, in the absence of zero-momentum modes, 
the classical vacuum at ^ = is isolated and the small volume behavior for the 
glueball masses is described by a perturbative series in g^{L), as opposed to 

Jence will be quite different, this 



Nevertheless, in large volumes 



g^/^{L) in absence of twist. The volume depcc 
is in particular true for electric flux energiest 
the results should not depend on boundary conditions. Therefore, comparing 
the different boundary conditions gives valuable information about the tran- 
sition to large volumes. In the hamiltonian formulation the twisted boundary 
conditions are most easily implemented in a gauge where the so-called twist 
matrices fij eSU(A^) are constant 

A(x-HLe(^)) = f7U(x)^-. (31) 



They satisfy 't Hooft's consistency condition^ which also gives the relation to 



the magnetic flux m S 'S?j. 



VL\.ri\VL^VLf^ — exp{2TTiekejrnj / N) . (32) 

These generate a so-called Heisenberg group (the group commutator is central, 
i.e. commutes with all group elements). The finite group theory allows one to 
construct in an elegant way the most general set. 17 ^ gS\J{N) for any given m, 
and its generalizations to higher dimensionsE3'E2l 

It seems that twisted boundary conditions spontaneously break the gauge 
invariance, due to the explicit choice of fife. However, this is of course similar 
to the case of periodic boundary conditions, which also represents a gauge 
choice in formulating the boundary conditions. Once one has specified the 
gauge choice for the boundary conditions, Gauss's law tells us that local gauge 
transformations have to satisfy 

h{x + Ln) = J|Ad"^(f7j)/i(x), Ad{n){h) = n^hn, (33) 



for n e TZj^, indicating the shifts over multiple periods. Note that in the adjoint 
representation the twist matrices commute. Thus, in a sense L is (typically) a 
1/A^ period (intimately related to the notion of color momentum, underlying 
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the principle of the Twisted-Eguchi-Kawai one-point lattice modelEa) How- 
ever, for finite size effects in large volumes, where the degrees of freedom that 
propagate "around" the boundary are colorless, this has no consequence (see 
Sec. 13). 

The candidate classical vacuum configuration satisfying the twisted bound- 
ary conditions is A = 0, of zero classical energy despite the presence of mag- 



netic fluxOO As we argued in Sec. |3.2| , the Polyakov loop (before taking the 
trace) , evaluated at the classical vacuum configuration should be invariant un- 
der gauge transformations. In the case of periodic boundary conditions, this 
implies it should be in the center of the gauge group, uniquely specified by 
Pj € Zff. To address the same question with twisted boundary conditions, the 
proper definition of the Polyakov loop has to be usedE3 



^.W^-tr 



tr (^P exp (^ ds Aj (x + se^^^ ) j n] 



(34) 



(path ordering from left to right) which indeed is invariant under the gauge 
transformations, Eq. (^3|). For ^ = this is even so before taking the trace, 
such that gauge invariance is indeed not spontaneously broken. Thus, Pj{A — 
0) = tr {il^j)/N and using Eq. ( p2[ ) one finds that Pj — 0, whenever nij ^ 
(mod N). This is closely related to invariance under constant gauge transfor- 
mations that are compatible with the allowed twisted gauge transformations, 

/ik(x + Ln) ==exp(27rik-n/7V) J|Ad"^ (f7j)/ik(x). (35) 

j 

As for the periodic case, k S Zf^ specifies the homotopy type of the gauge 
transformation. These gauge transformations multiply Pj with exjp{2TTikj /N). 
We note that A = is left unaffected by constant gauge transformations, 
/ik(x) = r2o(k), which from Eq. ( ^ have to satisfy 



17j(k)r2jr2o(k)f7^ = exp{2niki/N), (36) 

extending Eq. (|3^ ) to four dimensions. This equation is solvedJoLexample by 
rij, with kg — m,ieij£. In general solutions exist if and only iiEju k • m = 
(mod N). When k • m 7^ (mod N), the gauge transformation hy^{x) does 
not leave ^ = invariant, but maps to a vacuum state with fractional Chern- 
Simons number, see Eq. ( ^ ) (equal to i^(ft.k) as defined in Eq. (^5|)). It is 
separated from ^ = by a classical potential barrier related to the instaaton 
with fractional topological charge for twisted gauge fields on the torus to 
be discussed in more detail in Sec. |[ In Fig. |ll| we illustrate these features. 
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Since electric flux quantum numbers are associated with the representations 
of the Z]^ homotopy, this meana£3 that some of the electric flux states will 
have energies that do not vanish in perturbation theory, whereas the electric 
flux energies associated with tunneling through the classical barrier will be 
suppressed by exp(— STr^/A^g^). Both differ from the behavior we observed 
when m = 0. For example, for SU(2) and m = (0,0,1), with fii = iri, 
il2 — iT2 and il^ — 1, one finds ilo(l,0,0) = ±ir2, 17o(0, 1,0) = ±iri and 
no(l, 1, 0) = ±jr3 as the only non-trivial constant gauge transformations that 
leave A = unchanged. Therefore, energies of electric flux with ei and/or 62 
non-trivial are perturbatively lifted, whereas states that differ only by the 63 
quantum number are degenerate. 



L.V(A) 




Figure 11: Sketch of the SU(2) potential energy with twisted boundary conditions. Shown 
are two isolated vacua, separated by a barrier with height 34.148/ Lg^. These vacua are 
related by twisted gauge transformations /11/2 = hp with p - m = 1 (mod 2), having winding 
number !^(/ii/2) = 5- Gauge transformations = hi^ with k • m = (mod 2), leave the 
vacuum invariant. 

For SUI2I, the boundary conditions can be solved by the following Fourier 
expansionE3o 

A"(x) ^ ^ A'^(k-hr(''))exp(27ri(k-f r(")) -x/L), (37) 

with the color dependent "fractional" momentum 

= (eW X m)/2. (38) 

Momentum conservation ensures that gauge invariance is not broken by these 
color dependent momenta. Using the Coulomb gauge diAi — 0, one checks the 
classical vacuum A = is isolated with all fluctuations quadratic. To illustrate 
the perturbative analysis we consider the simpler case with m = (1,1,1), 
realized by flk = iTfc, because this choice of m does not break the invariance 
under the cubic group. Glueball states can thus be classified as for the periodic 
case. Generalizations to SU(3) (or arbitrary SU(iV) and magnetic flux) are easy 
to obtain. The allowed non-trivial constant gauge transformations are now as 
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follows: f7o(0,l,l) = ±iTi, f7o(l,0, 1) = ±iT2 and f7o(l, 1, 0) = iirg. The onc- 
particle state associated with the Fourier mode (p) , with creation operator 
6^(a,p, ±), has non-zero electric flux. They are further characterized by the 
momentum 27rp/-L, energy 27r|p|/L and the polarization ± of the gauge field, 
satisfying A^iji) • p = 0. To be precise, the electric flux vector belonging to 
this state is e = e^"^ This motivates interpretingo 27rr(°^/L as a Poynting 
vector. It also plays an interesting role in how the wave fpiictional behaves 
under translations over L in the three coordinate directionsLj For any 1^)^^ ^ 
with magnetic flux m and electric flux e 

|vI'(x + Ln))^^, =cxp(*n.Pi)|vI/(x))„^3, P = ^(e x m)/i. (39) 
The electric flux e = e'"^) is created with the gauge invariant Polyakov loop 



operator Pq, see Eq. (34). This contains the one-particle state h* {a,\),±)\Q) , 
the energy of one unit of electric flux is therefore given by the length of the 
Pointing vector, the minimal value 27r|p|/L can take, 

= 7r|e('') X m\/L = ttV2/L. (40) 

The energy of two units of electric flux (e.g. e = (0,1,1)) is perturbatively 
degenerate with this. At higher order one has to take into account that the 
two one-gluon transverse polarizations are now no longer degenerate. Properly 
creating electric flux with the gauge invariant Polyakov loop operator picks out 
the polarization in the direction of the loop; for the symmetric torus along e. 
This causes a jperturbative splitting between the energy of one and two units 
of electric fluxllj The energy ofpthxee units of electric flux, e = (1,1,1), is 
entirely due to instanton effectsp3'L3 see Fig. nl]. Therefore, in lowest order 
R2 = l + O and i?3 = -h C (exp(-47rVg^)) , quite distinct from what 
one finds with periodic boundary conditions in small and intermediate volumes. 

To find the mass gap in the zero electric fiux sector, one needs two-particle 
states, built from states with opposite (which for SU(2) is equivalent with 
identical) electric flux. They are of the form 

6t(a,p,«)6t(a,q,^)|0), (41) 

with a^(3 — ± the polarizations of the one-particle states. These states have 
total momentum Q and energy E satisfying 

Q = 27r(p + q)/L, - 27r(|p| + |q|)/i. (42) 

The minimal zero-momentum state gives the mass gap, in lowest order 

rugap = 27r|e(") x m|/L = iirV^/L, (43) 
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Table 1: One loop coefficientslljCa for SU(2) with m = (1, 1, 1). 
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which is twice the length of the Poynting vector. Counting the number of 
ways one can form these two-particle colorless zero-momentum states from the 
one-particle states, one finds a 24 fold degeneracy. This degeneracy will be 
lifted in. one loop order, arranged in irreducible representations r of |tlae cubic 
groupfa for which lattice discretization effects were reported as well|l3 but no 
details have been published. In the continuum, the O (g^) mass and energy 
shifts are parametrized by the constants 7^, 

Zr = rUrL = 2\/27r+^4S^' = E,^L = V27r+^47^ 3), (44) 

167r^ 167r^ 

which are listed in Table 0. 

In comparison to the case of periodic boundary conditions we note that 
the E^ tensor state is now heavier than the scalar At , but with the T2 in 
between. Also, is here a decreasing function of the volume, which has to 
turn around at some point, when the mass stabilizes and zq becomes linear in 
L. A clear finite volume artifact is also the near degeneracy of the oddball [A^ ) 
with the glueball {At). Both of these features we will also find for the sphere, 
see Sec. ^ There it will be demonstrated that taking the non-perturbative 
effects of instantons into account will lead to an appreciable splitting between 
the oddball and glueball. Also here, like for the case of periodic boundary 
conditions, we can estimate at which volume instantons become important, by 
equating the energy of the scalar glueball state mth the height of the barrier 
between two vacua, set by the sphaleron energy^ 34.148/L(7^(i). Again one 
finds the critical value of L to be of the order of 6 times the scalar glueball 
correlation length. First lattice Monte Carlo resiLl±s_,with twisted boundary 
conditioiis were obtained by Stephenson and TeperOS They find in very small 
volumeska (/3 = 4/gg = 4.7, on a 4^ x 96 lattice) that the Af , E^, and Tj+ 
glueball masses indeed all become degenerate and equal to 2Ee, with R2 = 1 
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and i?3 = 0. The A2 and states are appropriately heavier. Because the 
shifts in the masses are so small, a detailed comparison with the predictions 
for the O (g^) shifts is inconclusive. At larger volumes, Zq > 6, the Monte 
Carlo results show thatJiie differences between twisted and periodic boundary 
conditions disappearEZlEa 
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Figure 12: Numerical dataEd for the SU(2) electric flux energies, Ee^ = L'Eny/n, with 
twisted boundary conditions, m = (1, 1, 1), on lattices of size X Nt, with N = A,6,8 and 
resp. Aft = 128, 128, 64. The scale is set by L = 400Af exp(-3.38/3) fm and the curves show 
the perturbative predictions combined with a fit to the expected tunneling contributions. 

Also an extensive study was made of the electric flux energies 0@ Here 
lattices of size 8^ x 64 and N"^ x 128 with TV = 4, 6 were used, with /3 = i/g^ 
ranging from 2.25 to 2.6, reaching volumes between 0.3 and 1.3 fermi. In 
Fig. ^ E„ = £'e„ / (Ly/n) is plottecLas a function of L. This definition is 
such that if, as predicted by 't Hooftp i?„ ^Jn in the infinite volume limit, 
S„ gives the infinite volume string tension K independent of n. The curves 
test the small volume expansion. J^hey contain the perturbative contribution, 
corrected for the lattice artifactsO (the data is not accurate enough to test 
the O (go) corrections), together with const. xg^^exp(—5'/gn); the expected 



shift due to tunneling (cmp. Fig. |TT| and Sec. 6^) mediated by the charge 
i instanton with classical action S {An^ in the continuum, see Sec. ^). The 
constant is fitted separately for = 4 and 6 (A^ = 8 fits not shown). One 
finds fair agreement with the predicted behavior, with confirmation that the 
transition to the large volume starts around 0.75 fermi. 
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3.5 Super symmetry and the Witten Index 

The case of supersymmetric Yang-Mills theopes in a finite volume has been 
considered in the context of the Witten indexlld in some detail. The torus ge- 
ometry is crucial to preserve the gauge invariance. The Witten index involves 
counting the number of quantum states (fermionic states with a negative sign). 
At non-zero energy, because of the supersymmetry, this number cancels be- 
tween the bosonic and fermionic states, such that the counting can be reduced 
to the vacuum sector. Here the number can be non-zero, because the super- 
symmetry generator can annihilate vacuum states. A zero Witten index is a 
sign of spontaneous breakdown of supersymmetry, where the vacuum energy 
is non-zero, explaining the physical significance of this index. A zero vacuum 
energy is a direct consequence of unbroken supersymmetry, where the hamil- 
tonian is given hy H = ^{QjQ''} with Q, the supersymmetry generators, 
that annihilate the vacuum state. 

In perturbation theory bosonic loops are typically cancelled by fermionic 
loops, e.g. in the vacuum energy. Applied to the problem of non-abelian 
gauge theories in a finite volume this implies that the vacuum-valley effective 
potential vanishes. The cancellation is caused by the contribution from the 
gluino fluctuations, which are the superpartners of the gluons. They are Weyl 
fermions in the adjoint representation of the gauge group, denoted by X^, 
with a a two-component spinor index. This means that the wave function is 
no longer localized to A = or any of the other orbifold singularities of the 
vacuum valley (the moduli space of flat connections) . It has been shown in the 
context of the supermembrane, when taking the supersymmetric Yang-Mills 
hamiltonian restricted to the zero-momentum modes only, that the spectrum is 
continuous, down to zero-energy. One can construct trial wave function&.with 
a support arbitrarily far from A — that nevertheless have finite energyO The 
compactness of the vacuum valley is crucial to obtain a discrete spectrum. 

The counting of quantum vacuum states was based on the assumption 
that for all gauge groups the moduli space of flat connections is given in terms 
of the Cartan subalgebra, as we discussed for SU(A^). The gluonic part of 
the groundstate wave function |0) is assumed to be constant over the vacuum 
valley. In the reduction to the vacuum valley there are r gluinos (associated 
with the generators in the Cartan subalgebra), each with two helicities. They 
are constant and carry no energy, which is the source of the vacuum degeneracy. 
These gluinos have to be combined in Weyl invariant combinations, respecting 
Fermi-Dirac statistics. There are r independent invariants, made from 

U = Sabe'^^XlXl (45) 

and its powers. So one has C/"|0), n = 0, 1, • • • , r, as the r + 1 bosonic vacuum 
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states, with no invariant fermionic vacuum states. Thus, one finds an index 
equal to the rank of gauge group plus one, r + 1. 

Because of possible problems with the adiabatic approximation near the 
orbifold singularities, as we encountered in the previous^sections, Witten con- 
sidered the alternative of twisted boundary conditionsllEl For SU(A^) the same 
result, r + 1 = N, for the index follows. Other groups in general do not admit 
the type of twisted boundary conditions that completely remove the continuous 
vacuum degeneracy with its orbifold singularities. So it was natural to attempt 
to find the exact zero-energy ground state for the supersjymmetric generaliza- 
tion of Liischer's zero-momentum effective hamiltonianp^ or even ignocing the 
time dependence, studying the path integral in the ultra local limitB3 None 
of these studies took the compactness of thej-uacuum valley into account and 
therefore fail to address the proper situationu3 

The problem of the adiabatic approximation remained an urgent one be- 
cause of a discrepancy between the finite volume calculation of the Witten 
index, and the one based on the infinite vqliirnft-determination of the gluino 
condensate through instanton contributionsE30'Ej One relies on the fact that 
an index can not change under smooth perturbations, like increasing the vol- 
ume. Also the infinite volume calculation has its problems. It uses the semi- 
classical approximation for a strongly interacting theory. This could, however, 
be circumvented by first adding matter fields to introduce an external mass 
scale to control the instanton calculation and then rely on the index being 
constant under a smooth deformation (through holomorphy), that decouples 
the extra matter sectorEZl This resulted nevertheless in a discrepancy, ■y/5/4 for 
SU(2), between the so-called strong and weak coupling calculations. Both cal- 
culations rely on the cluster decomposition property, since the instantons have 
more than two gluino zero modes, which seems to make the condensate (AA) 
vanish. The instanton calculation instead computes the appropriate power of 
the gluino condensate, ((AA)''), that saturates the 2h gluino zero modes, where 
h is the so-called dual Coxeter number of the gauge group, h = N for SU(iV). 
It is this power that gives the number of vacuum states, 

(AA) EE e^""/'' (I ((AA)'') I) n= l,2,---,/i. (46) 
These arguments seepLreasonable, but are not rigorous^ See for further details 
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been made of the constituent nature 
in the context of a Kaluza-Klein re- 



a review by ShifmanO Recently, 
of periodic instantons (or calorons) 
duction with periodic gluinos, as opposed to a high temperature reduction with 
anti-periodic gluinos, which would break the supersymmetry. The constituent 
monopoles have exactly two zero-modes and saturate the condensate, (AA). 
The strong coupling calculation now agrees with the weak coupling result 
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The period can of course be used to control the couphng constant, but it is 
assumed the index does not change, going from a small to a large period. 

The mismatch in the Witten index between small and infinite volumes 
occurs for SO{N > 6) and the exceptional groups. There has, however, been 
a recent revision in counting the number of vacuum states pin, a finite volume. 



In a study of D-brane orientifolds in string theory, Wittenll23 constructed for 
S0(7) an extra disconnected component on the moduli space of fiat gauge con- 
nections, which can be embedded easily in SO{N > 7). For S0(7) and S0(8) 
this gives an isolated component of the moduli space, contributing only one 
extra vacuum state. For SO(A^ > 8) the extra component in the moduli space 
behaves like the trivial component for S0(iV-7). Adding r + 1 coming from the 
SO(A^) and SO(A'^— 7) moduli space components gives the dual Coxeter num- 
ber of SO(A^), thereby giving the same number of vacuum states as obtained 
in the infinite volume. 

Witten's construction based on orientifolds does not work for the excep- 
tional groups. This natuiplly led to a derivation of the extra vacuum states in 
a field theoretic contexttSil trivially extended to the exceptional group G2, 
as a subgroup of SO (7). Three different groups have independently man- 



aged to sohie. the problem for other exceptional groups with periodic 



conditionai-J'tii^l and for any group with twisted boundary conditions 



iundary 



!£Zl As we 



remarked before, twisted boundary conditions usually do not remove all the 
vacuum degeneracies, but it is important that the number of vacuum states is 
independent of the twist for all gauge groups that have a non-trivial center. 
The origin oLthe extra moduli space components is actually not too hard to 



understandE£3 Large gauge groups can have subgroups that are products of 
unitary groups, which each would allow for twisted boundary conditions. By 
choosing twists from all subgroups to cancel one obtains periodic flat connec- 
tions that can not be deformed to the Cartan subalgebra, which supports the 
trivial componenLof fiat connections. Of course one need not cancel these 
twists completelyE£] Needless to say, the group theory involved to sort out all 
the constraints and count the number of vacuum components is rather involved. 

Supersymmetry does not play a role in establishing the existence of these 
extra vacuum components. Supersymmetry is, however, crucial for these ex- 
tra components to lead to extra quantum vacua. As soon as the perturbative 
quantum fluctuations in the vacuum energy do not cancel, this will in general 
be different for different vacuum components. In a small volume, i.e. at weak 
coupling, the wave functional will localize around the one with the lowest vac- 
uum energy. Within such a connected component it will localize around the 
minimum of the effective potential, as we discussed in the previous sections. 
It is likely, since the trivial vacuum component is the widest, that this is the 
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one where the wave functional locahzes. But interestingly, one now has po- 
tential energy barriers between vacuum components that are not related by 
a homotopically non-trivial gauge transformation (since the different vacuum 
components are not isomorphic) . Still these vacua can be characterized by frac- 
tional Chern-Simons numbers and tunneling between them would be dsacxibed 



by new types of instanton solutions with fractional topological chargeliiU'tlij It 



considerably adds to the richness of non-abelian gauge theories. 

Although these new results for counting the number of vacuum states in a 
finite volume remove the urgency of addressing the problem with the adiabatic 
approximation, it does remain a wnrp point in the finite volume analysis, as 



also stressed recently by WittenE23 One immediate problem we encounter is 
that Eq. (^ ) seems to imply that the vacuum-valley wave function has to 
vanish at the orbifold singularities, which seems inconsistent with it being 
constant sufficiently far from the orbifold singularities. However, one should 
take into account the behavior under Weyl reflections of the occupied negative 
energy gluino states in the Dirac sea near the orbifold singularities in the 
zero-momentum hamiltonian. Attempting to incorporate this in the formalism 
developed in the previous sections, we run into the problem that spin i fields do 
not decompose in three components, which can be associated with each of the 
coordinate directions. In the bosonic sector we could conveniently ignore the 
transversality, with gauge invariance restored at the end by restricting to the 
invariant wave functions. Thereby each of the coordinate directions separately 
allowed a polar decomposition, not compatible with the nature of the gluino 
fields. For SU(2) there is, however, an elegant polar decomposition Jising the 



spherical and gauge symmetry of the zero-momentum hamiltonianPi^2l which 
even holds for Hcs in Eq. (^) to O (5^/^). In the supersymmetric case one 
would require the wave function to become constant after reduction to the 
vacuum valley, sufficiently far away from A = Q, to match to the expected 
behavior away from each of the orbifold sinsuJarities. Such a matching can in 
principle be controlled sufficiently rigorously^ as was done in the semiclassical 
calculation of the energy of electric fluxplEJ but here we do not have the benefil 
of a well localized zero-momentum wave function atlhe orbifold singularitiesp3 
and judging the incomplete results in the literatureEj this seems not the way to 
go. It is in the light of the robustness of supersymmetry somewhat surprising 
(and frustrating) this problem remains so technically demanding. 



4 Instantons and Sphalerons on the Torus 

Instantons are associated with the tunneling through barriers that separate 
different vacuum components. We are interested in these barriers to establish 
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the directions in field space in which the spreading of the wave functional has 
to be taken into account, when energies of the low- lying states no longer are 
well below the barriers. The barrier is a saddle point of the energy functional 
with one unstable mode, a so-called sphaleronEil They exist here by virtue 
of tha^finite volume. Remarkably, on any compact four manifold instantons 
existp3S still with an arbitrary scale parameter, only limited by the size of 
the manifold. The finite volume sphaleron is typically associated with the 
instanton that achieves this maximal scale. 

For low topological charge, sometimes smooth instantons can not be proven 
to existE3 An example of non-existence actually occurso on T'*. In general 
one takes a small localized instanton from that is matched smoothly to the 
fiat connection. When the moduli space of fiat connections has a continuous 
degeneracy, as is the case for the torus, there can be obstructions against this 
procedure. One can get arbitrarily close to satisfying the self-duality equations, 
Ffii, = iF/^ii, = ^^e^vapF'^^ , that saturate the bound for the euclidean action 

S=-i J tr{Fp = -i J tr{F^,±F^,f±i J tr(F^,F^'') <87r2|z/|, (47) 

with v the topological charge. However, an exact self-dual solution is only 
achieved in the limit where the scale parameter is forced to zero or the volume 
is taken toiiifinity. For twisted boundary conditionsLI, existence of charge one 
instantonsO can be understood from the fact that twist removes the continuous 
degeneracy in the moduli space of fiat connections. 

From the hamiltonian point of view we are interested in instantons on 
T'^ X M, where time is not constrained to be finite. This can be seen as a limit 
of T'' where one of the periods, called T, is taken to infinity. In the Aq — 
gauge the field is periodic up to a gauge transformation h in the time direction, 
A(x,r) = [h]A{x,0), and 
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•T 

tr {F^.F'^'') = / CS{A) = CS{[h]A) - CS{A) = v{h), (48) 



with v{h) defined in Eq. (|T|). Using that tr {F^^F'"') = d^,K^{A), we find for 
the Chern-Simons functional 

CS{A) = J d^xKoiA) = Jd^xtY {A,{d,Ak + lA,Ak))e,,k- 

(49) 

When has periodic boundary conditions there is an ambiguity which point 
in the vacuum valley to tunnel from and to. These need not be the same points, 
as is already evident with tunneling related by non-trivial twisted gauge trans- 



fomations (see Eq. (20) and (35)), under which the Polyakov loop is periodic 
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in time up to a non-trivial Zj\j_factor. In that case charge one instanton so- 
lutions can be shown to exist &3 The obstruction for the periodic four-torus 
is reflected in the fact that the instanton does not want to tunnel between 
the same points in the vacuum valley. This is illustrated in Fig. K^, obtained 
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Figure 13; Numerical result£3 obtained with cooling on A'^'^ X 3N lattices (scaled to L = 1). 
Left is with k = 0, A'^ = 8 and right with k = (1,1,1), 7V = 7 and 8 (k indicating the twist in 
time). The top figures show electric (££;(t) triangles) and magnetic (^^(t) squares) energies 
(g = I), the inset shows the tails at an enlarged scale. In the lower figures we plot C2(t) 
through two distinct (on the right one) point (s). 



by a numerical procedure called coolingE_l that minimizes the lattice action 
in a giveiutopological sector, carefully dealing with the lattice discretization 
problems^ With periodic boundary conditions the configuration first relaxes 
to a solution which is not self-dual, although it has reached the vacuum valleys 
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Table 2; The lowest eigenvaluesc2l for the Hessian of the energy functional evaluated at the 
twisted and periodic sphaleron on a A'^'^ lattice for = 4, 6 and 8. 



m=( 1,1,1) sphaleron 


m= 


=0 sphaleron 


N=4 


N=6 


N=8 


N=4 


N=6 


N=8 


-11.1280 


-10.9003 


-10.8174 


-13.1960 


-13.3935 


-13.4432 


0.0074 


0.0001 


0.0000 


0.0073 


0.0006 


0.0000 


0.0091 


0.0001 


0.0000 


0.0124 


0.0010 


0.0000 


0.0101 


0.0001 


0.0000 


0.0139 


0.0011 


0.0000 


13.1667 


13.3371 


13.3934 


1.7041 


1.4151 


1.3009 


13.1667 


13.3371 


13.3934 


1.7070 


1.4303 


1.3012 


13.1667 


13.3372 


13.3934 


7.9677 


8.0984 


8.1260 


16.0475 


16.8663 


17.1597 


7.9682 


8.1008 


8.1261 


16.0478 


16.8663 


17.1597 


7.9684 


8.1279 


8.1266 



at either end in euclidean time, judging the zero magnetic energy. However, 
it moves in accordance with the equations of motion along this vacuum valley 
with constant speed (i.e. electric field) to match the forced periodicity, as is 
also evident from the linear (space independent) slope for C2(t). When one 
would continue to lower the action, in order to reach a self-dual solution, the 
instanton is forced to shrink to zero size. With a twist in the time direction 
the configuration immediately relaxes to a self-dual one. 

Shown on the right in Fig. |l^ is an instanton that is close to being maximal 
in its size, thereby going through the lowest barrier, to maintain the same 
action. ThiaJowest barrier is assumed to be a finite volume sphaleron. This 
was verifiedp3 by explicitly looking for saddle-point solutions on T"^, minimizing 
the square of equations of motion using a similar numerical procedure of cooling 
(cmp. Fig. ^|). This saddle point is not invariant under translations, giving 
rise to zero-modes, which in addition to the unstable (i.e. tunneling direction) 
will become the degrees of freedom that need to be treated non-pcrturbatively, 
by imposing the proper boundary conditions in field space. In addition it was 
found that there are two directions not associated with zero-modes, in which 
the saddle point is nearly flat (see Table ^). Such modes would have to be 
treated non-perturbatively as well. Without an analytic solution, it is hard to 
judge if such an analysis is feasible. This was one of the motivations to look at 
the situation where the torus geometry is replaced by that of a sphere. In this 
case the finite volume sphaleron has a constant energy density and no (almost) 
flat directions, see SecH 

We already mentioned in the previous section that with non-zero magnetic 
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flux, there exist instantop_solutions with fractional topological charge. When 
first proposed by 't Hooftl3 this was met with some disbelieve. The mathemat- 
ical classification of so-called fiber bundles, to describe gauge fields on compact 
manifolds seemed to require integer topological charge, specified by the second 
Chern number (see Eq. (^J)). But due to the twist, the fiber bundle has to 
be considered real and the appropriate classification is by the first Pontryagin 
class. With the proper normalizations, this makea^the would-be fractional sec- 
ond Chern class an integer first Pontryagin classEj It is the difference between 
a unitary and orthogonal fiber bundle, e.g. SU(2) versus SO (3). This makes 
also precise the intimate connection between the topological charge and twist, 



as first discovered by 't HooftO As we have seen in Sees. and 3.4 (see 
Eq. (|3^), which also applies for m = 0, in which case flj = 1), there are ho- 
motopically non-trivial gauge transformations /ik, "periodic" up to an element 
of the center of the gauge group, that do not affect the boundary conditions 
of the gauge fields. In particular h^^dh^^ has the right periodicity, and v{h-k) 
in Eq. ( p^ ) remains well-defined (in the gauge where the twist matrices Vtj are 
constant and Aq = 0). When m = we have seen that we can always choose 
/ik abelian, and i^(ft.k) — 0. In this case there is no interplay between the twist 
and the topological charge which is here integer. When m is non-trivial we 
stated in Sec. |3.4| that when k • m = (mod N), /ik can be chosen constant, 
= f^o(k), as can be checked for the simple case of SU(2) by inspection, but 
can be proven from first principlesEZI However, when k • m 7^ mod N it may 
be that i^(ft.k) 7^ 0. The additive nature of the topological invariant ^{h), 

v{h'h) = v{h') + v{h), (50) 

implies that v{h^) — Nv{h\f^). This has to be an integer, since is "periodic" 
and the winding number of a periodic function /i : ^ G is an integer. Also 
one can duplicate , under which the topological charge is additive, as many 
times in the space directions as is necessary to get a newJxjrus for which m = 0, 
and hence with integer topological charge. ThereforepI v{h) is proportional 
to and linear in k and m, such thaLj^ + m • 'k/N G TL, as was proven 
rigorously in the context of fiber bundlesLiJ Conversely, using this result one 
easily sees that the existence of constant twist matrices in four dimensions, i.e. 
solutions to both Eq. ( |3^ ) and Eq. (|36|), require k • m = mod N. This is 
so because ^ = 0, with = 0, satisfies the boundary conditions. For further 
details on these twist-related issues one can also consult a recent reviewE3 

The above implies that the lowest non-trivial topological charge that can 
be reached is 1/iV. With Eq. (^7|), the associated euclidean action of this in- 
stanton solution is found to be Stt^ /N, which can indeed be attained as one 
can prove existence of self-dual solutions in the appropriate sectorEa Like for 
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topological charge one, despite much effort up to date nQ_analytic solutions 
are known. What can be deduced from general principlea23 is that the solu- 
tions with charge \/N have only the four translational degrees of freedom. In 
particular they have no free scale parameter, their size is fixed with respect 
to the torus. JXhis-fixed scale is a great benefit in studying these instantons 
on the latticeoEjO For this fractionally charged instanton the finite volume 
sphaleron was constructed as well, by finding for SU(2) with twist m = (1, 1, 1) 
on the appropriate saddle point. Like for m = it is not invariant under 
translations, giving three flat directions at the saddle point, but in this case 
no other near zero-modes appearpi see Table. | 



S 8 




Figure 14: Numerical results^j for S{t), the square of the equations of motion, along the 
tunneling path for the SU(2), charge i twisted instanton solution for m = k = (1,1,1), 
evaluated on a 8^ X 24 (triangles) and a 12^ X 36 (squares) lattice (scaled to L = 1). The 
dip at t = occurs at the finite volume sphaleron. 



In Fig. |l^ we illustrate the fact that along the tunneling path there is a 
saddle point, by plotting the square of the equations of motion. 



S{t) 



1 

32 



T3 



<fx tr {d,F,j (x) + [A,(x), F,, (x)])^ 



(51) 



as a function of euclidean time. The fact that S does not exactly vanish in the 
middle, at the top of the barrier, is due to a limited resolution on the finite 
lattices employed. 
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T- duality for Instantons on the Torus 

Although no explicit solutions of the basic instantons of charge 1/iV and charge 
1 are known, progress has been made in better understanding instantons on the 
torus. We will first discuss the case of T'*, before addressing x TR. For 
the only explicitly kno'^yn.soluti.ons are the ones constructed by 't Hooft with 
constant field strengthEi!Jlliillii3 These are essentially abelian and are self-dual 



only when the periods of the torus match the field strength (which is quantized 
due to quantization of flux). Recently these solutions to the self-duality equa- 
tions were_&pllowed in perturbation theory by deforming away from th£_special 



geometryfij using as a starting point the exact fluctuation spectrumMll in the 



background of the non-selfdual constant curvature solutions. 

The essential ingredient for gaining fm±h£r understanding has been the 

Nahm introduced it to find 



so-called Nahm or Mukai transformationtl 111 



magnetic monopole soliitipns, getting inspiration from the algebraic Atiyah- 
Drinfeld-Hitchin-Manir£i3 (ADHM) construction of instantons in IR^. The 
formalism can as a mailer of fact be used to provide a simple proof of the 
ADHM construction E3Ei3 It will be convenient in the following to see as 
IR^/S, where S is a four dimensional lattice spanned by the four periods a^'^\ 
and to introduce the connection one- form lu{x) = A'^{x)dx^, which is invariant 
up to a gauge transformation under translation over a period. These gauge 
transformations, called cocycles, satisfy cocycle conditionsQ (twisted boundary 
conditions will be discussed later) , to assure the gauge invariant quantities ace 
periodic, and such that one has an appropriate fiber bundle over the torusp3 

uj{x + a) = h^^{x){uj{x) + d)ha{x), (52) 
ha+bix) = hb{x)ha{x + b) = ha{x)hb{x + a), a,beE. 

It will be convenient to take V{N) as gauge group, allowing for a U(l) factor. 
The associated vector bundle in the fundamental representation of U(-/V) is 
denoted by E. The curvature of this bundle E, given by 

F da; + UJ ALU = ^Ff'^dx,, Adx^, (53) 

leads to a topological charge ly — ^^^{c2{E) — ^ci{E) A ci{E)) as an integral 
over Chern classes (cmp. Eq. ( p8| ) and note that ci{E) = for SU(A^) because 
of the vanishing trace), where 

c.(i.) =tr(^), c.(i.)^.tr(Z.,Z.). (54) 

The starting point of the Nahm transformation is that gauge fields with 
topological charge i^, have v positive chirality zero-modes for the massless Dirac 
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equation^ which is most conveniently written in the Weyl representation. For 
this we introduce the unit quaternions cr^ and their conjugates = cr^ 

(Tp (l2,irj), 5-^ = (I2, -iTj). (55) 

They satisfy the multiplication rules 

O'p^!^ = V'jl,y<^a, ^^CTi' = ^^lyCTa, (56) 



where we used the 't Hooft 77 symbolsp3 generalized slightly to include jyJJ^ = 
fj^^ = 5^y, also useful in Sec. ^. The Weyl-Dirac operator is given by 

D = a,,D^{A) - (7^(9,, + A,,). (57) 

Hence, in the background of a charge v gauge field there are v independent 
solutions to = 0. For 'i>{x) to be defined as a two-spinor on the topji 
one requires '^{x + a) = haix)"^ (x) . Strictly speaking, the index theorem 
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only states that the number of zero-modes of positive chirality minus those 
of negative chirality (for which D^vj/ = 0) is equal to the topological charge, 
but it will be assumed that there are no negative chirality zero-modes. For 
self-dual gauge fields we will discuss this assumption later. One nnw adds a 
spectral parameter z £ in the form of a flat abelian connectionliiJ 

uJz — Lu + 2'Kiz'^dx^, (58) 

(the unit generator of U(l) is implicit) which leaves the curvature unchanged, 
Fz — F . Hence there is a smooth family of u normalized fcrmionic zero-modes 

D,^^\x) = a^(a^-f v4^ + 27rzz^)*W(a;) = 0, / d^x ^^\x)^^^^\x) = 5,,. 

(59) 

From this family one constructs the connection Qj = Ap(z)(izp by 



i;^(z)= / d^x^l{x)^^^^p{x). (60) 

Using that i:)^+a(e"^''" '^^'i''(a;)) = and the fact that {^'i+a(a;)} forms a 
complete orthogonal set of solutions, we find that 

*i+a(^)=e~'""-'*i''n^)^rW, (61) 

with ha{z) a unitary v v matrix, which defines the cocycle for Q) 

Cj{z + d) ^ hl^{z){u}{z) + d)ha{z), d = dz^dz^, (62) 
45 



as a connection on the dual torus = M^/S, where 

S = {a e IR^'I < d,a >e Va e S} (63) 

is the lattice dual to S. For example, if 5 is generated by a^'^-' = Le'^^\ giving 
a hypercube with sides of length L, the dual lattice S is generated by d^'^^ = 
7j-ig(M)^ a hypercube with sides of length \/L. So the Nahm transformation 
is a T-duality. 

In the absence of negative chirality zero-modes (cokerDz = kaL,p|_= 0) 



Ez = kerZ?2 depends smoothly on z. It was realized by Braaml^lii3 that 
one can use the Atiyah-Singer family index theoremlli3 to compute the Chern 
character of the bundle E with connection uj, as the integral over T'* of the 
Chern character of E ® V, which is the vector bundle relevant for Uz, as 
a vector bundle over x T^. Here V is the so-callecO Poincare line bundle 
with the connection i'Kiz^^dx^. It was crucial this line bundle had no curvature 
as a bundle over T'', but extending it to T'^ x T^, this is no longer the case 
and Fjy — iiridx^ A dz^, with a first Chern class, Ci{V) = dx^ A dz^^ (for 
a line bundle all other Chern classes vanish). The Chern character satisfies 
ch{E (X) £") = ch{E) A ch{E') and is a formal polynomial in the Chern classes, 

ch{E)^co{E) + ci{E) + ^ci{E)/\ci{E)-C2{E), ch{V) = exp{ci{V)). (64) 

Perhaps somewhat confusingly cq{E) = rk{E) is called the rank of the vector 
bundle E, here the dimension of the fundamental representation (so rk{E) — N 
for a V{N) fiber bundle). It now follows that 

ch(E) = / ch(E) Ach(7'). (65) 

The integral over only gives a non-zero result if the combined form is of 
degree four in dx^, i.e. a top- form. This means that 

• The zero- form co{E) = N, becomes a volume form on T^, whose integral 
over gives the topological charge of the dual gauge field. 

• The two- form ci{E) whose cohomology class can be specified by the 
a&eZian fiuxes ci{E) = in^''da;^A(ia;,y, goes over in ci{E) — ^n^'^dZfj^Adz^. 

• The four-form whose integral is the topological charge v of the original 
gauge field becomes a zero- form, cq{E) — v. 

Thus we see that the Nahm transformation interchanges the rank N of the 
bundle with the topological charge and maps the first Chern class to its 
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dual. In particular, an S\J{N) bundle has vanishing first Chern class and 
is thus mapped under the Nahm transformation to U(z^) with vanishing first 
Chern class. It is well-known that in the latter case the bundle can be gauged 
to an SU(i/) bundle (its U(l) part is trivial and can be gauged away). 

The family index theorem only provides topological information on the 
dual gauge field. We will demonstrate the wondrous resultEJ, that if w is a 
self-dual connection, than also lu is self-dual. A crucial ingredient is formed by 

D,Dl = -DliA,) - (66) 

using cr^CTi, — 6^1, -t-o-j^a^]. Since cr[^a^] = CiTy^^ is an anti-selfdual tensoi^we 
see that D^D\ — —D'^^{Az). We first deal as promised with the condition that 
there are no opposite chirality zero-modes. We have seen that if oj is self-dual, 
and Dl-^ = Q that Dl^^A^)-^ = 0^-^ = 0, but this would imply that the 
connection oj has a flat factorJS meaning that the bundle of rank N splits in 
the direct sum of a bundle of rank TV — 1 and a flat line bundle, E ^ E' ® L^, 
where has the connection Such a flat factor would exist for any 

z (simply shifting ^ proportionally), and so a sufficient condition to impose 
the absence of opposite chirality zero-modes is to require oj to be without flat 
factors (WFF, sometimes also called 1-irreducible). 

A direct corollary is now that regular SU(iV) charge one instantons cannot 
exist on T^. Suppose they would exist. The Nahm transformation gives rise to 
a U(l) bundle of charge TV, which is impossible as the first Chern class vanishes 
and U(l) bundles have always vanishing second Chern class. Suppose flat 
factors prevent us from making the argument. In that case we can remove the 
flat factor and go down in rank by one (which keeps the first Chern class zero). 
This can be repeated until we are left with an SU(iV — k) charge one instanton 
without flat factors, or completely reduce to flat factors which, however, can 
not support topological charge. 

Assuming no flat factors for the self-dual connection w, the kernel of D^D\ 
is trivial and the well-defined Greens function Gz — {D^D\)~^ commutes with 
the quaternions cr^. It is now straightforward to show that lo is self-dual as 
well 

F'^{z) = duj'^ + uj'^ A uj*'^ dl-W |1 - Pz\Mp >, (67) 
where Pz is the projection on kerZ?^, 

Pz = Y^ l^'i''' >< = 1 - DlGzD,. (68) 
fe=i 



''Note tlmjt 't Hooft had introduced X4, instead of xq, as the imaginary time. Using his 
symbolsE3 replacing the index 4 by 0, flips the duahties around, since eoi23 = —£4123- 
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The equahty to the right-hand side is triviaUy true for a zero-mode, whereas 
the orthogonal complement of kerl^^ is described by the image of D^, on 



which both sides act as the identity. Next we use D^cM! 



to find 



(69) 



Self-duality immediately follows from the fact that CT^^cr^] — Oiff^^ is self-dual. 

For T^, applying the Nahm transformation the second time, it can be 
shown that lo is retrieved identically, and the explicit form of ^x{z) in terms 
of ^z{x) allows one to show that metric and hyperKahler structures of the mod- 
uli spaces are preserved under the Nahm transformation. In other words, if we 
denote by M.N,iy the moduli space of SU(iV) charge v instantons, the Nahm 
transformation induces a map between moduli spaces, N : Mn,i^ ~> ■A^j^.a', 
which is an involution that preserves the natural metric and hyperKahler struc- 
ture of the moduli spaccEJ The dimension of the moduli space, ANv, is indeed 
symmetric under interchanging v and N. 
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Figure 15: Shown is a cross section through the 1-2 plane of the various unit cells that appear 
in the Nahm transformation for SU(2) with twisted boundary conditions ni2 = — n2i = 1 
on i?*/S. We have chosen here Ha^^^^H = Ha^'^'H = ^V2. Duplicated unit cells without twist 
are S' and S", their union is denoted by So- Its dual, Ho, gives the torus on which the 
Nahm transformed gauge field lives, with the same gauge group and twist as for the gauge 
field we started with. The unit cell Ho is also obtained by intersecting H' and H", the duals 
of the duplicated unit cells. When a symbol overlaps with different cells it belongs to all of 
those cells as a whole. 



The Nahm transformation on can be defined to innhide twisted bound- 
ary conditions, as Gonzalez- Arroyo has shown recentlyt£3 The basic idea is 
simple: One duplicates the unit cell, also denoted by S, in the various direc- 
tions as many times as is necessary to remove the twist. Let us assume for 
SU(2) ni2 = —n2i = ms = 1 and the rest of the twist factors trivial. We 
can duplicate S either in the 1- or in the 2-direction in order to remove the 
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twist, giving ^' and In both cases we can apply the Nahm to the enlarged 



unit cells S' or S", see Fig. 15, In each case its dual S' or S" is shorter by a 
factor 2 in the direction we originally duplicated. These two choices of unit cell 
intersect in a cell we call Sq that is twice shorter in each of the two directions. 
Due to the intersection it contains all the gauge invariant information of Cj' and 
ill" . The gauge field therefore has half-periods associated with twisted bound- 
ary conditions on H'^/Sg- In this particular example the resulting twist is the 
same as we started with. The argument can be geiienalized, either along the 



above lines, or using so-called flavor multiplication^£3t£il to arbitrary SU(iV) 
and twist. The example shown in Fig. |l^ (with identical behavior in the 0-3 
plane), maps the charge i instanton to the same solution. Another example 
of self-similarity occurs for SUC21 in the charge 2 sector. When applied to 



the special solutions of 't Hooftliil] this can be used as an example where the 
Nahm transformation can be worked through analytica|lly. The resulting Nahm 
transformed connection is again of the special formlJ 



To have the Nahm transformation help us finding explicit instanton solu- 
tions, simplifications have to be arranged for. These happen to be perfectly 
geared to finding charge one instanton solution on x M. For charge one the 
Nahm transformed gauge field is abelian and with one period infinite, the dual 
period collapses to zero. The decompactification limit T — > oo is associated 
with dimensional reduction. This reaches a dramatic height for the case that 
all periods are send to infinity, relevant for instantons on IEI^ The dual is 
now a single point, explaining why the ADHM constructiorilij is in essence 
algebraic. However, under the decompactification limit, a partial integration 
is required in going from Eq. ( |67| ) to Eq. (^) , 

+87r2^^, <vI/W|afcG.|*i^'^ >. (70) 



lie 



These boundary terms destroy self-duality of F, which can be repairedti 
for instantons on IR'* and for calorons (periodic instantons) on M'^ x Si . Except 
for M^, for which there is no freedom, the location of the singularities are fixed 
by the asymptotically fiat connections at t — > ±cx3, which are required to occur 
to keep the action finite. The Weyl-Dirac hamiltonian in the background of a 
fiat connection has generically a mass gap. The mass gap vanishes, however, 
in particular when the flat connection is pure gauge. The flat connections are 
characterized by — 27ri diag(w(^^ + z, • • • , vif^^^ -I- z), and for the mass 
gap to vanish it is sufficient that \[w'^^^ + z|| = for any of the N values 
of k. With two asymptotically flat connections, their are 2N corresponding 
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values of z (which are defined modulo S) where the mass gap vanishes. Only 
here the boundary terms in Eq. ( |70| ) can be non- vanishing. A simple steepest 
descent analysis in their neighborhood shows that they act as point sources 
with abelian electric and magnetic charges ±7r, enforced by charge quantization 
(positive charges associated with the flat connection at < = oo). This ensumt 



the magnetic sources are Dirac monopoles with unobservable Dirac strings 
Outside of the singularities the abelian field is self-dual and, noting the zo 
independence, Bj{z) = Ej{z) = —dAo{z)/dzj, where 

. N 

For the usual normalization of abelian fields one multiplies this with —i. The 
sum over the periods S on T^, although formally divergent ,_caji be resummed 



in a rapidly converging series and one has an exact resultErS for Ao{z) (see 



Fig. |T|). One may argue from this that indeed time-periodic instanton so- 




Figure 16: Comparisont£j of in the plane Z3 = 0.5 for the Nahm transformation of an 
SU(2) instanton on an 8^ X 40 lattice (scaled to L = 1) and twist k = (1, 1, 1) in the time 
direction (left), with the analytic result on X IR for P~ = (0.86, 0.76, 0.08) = -P+ (right). 

lutions do not exist, as for w^-* = w^-* we clearly find Aq{z) ~ 0, which 
under the Nahm transformation leads to a trivial connection on x K. With 
twisted boundary conditions in the time direction one has w^'' = w'i''' + h/N. 
Although it is not straightforward to solve the Weyl-Dirac equation on in 
the background of an abelian self-dual gauge field generated by point charges, 
one may hope explicit solutions of charge one instantons on the torus can be 
found this way. 
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The situation with the twisted instantons of fractional charge is, perhaps 
somewhat paradoxicaUy, a bit more comphcated. However, much understand- 
ing has been gained through the numerical implmientation of the Nahm trans- 
formation on a latticetHil It has been foundlH^ that, in staying away from 
the exact decompactification limit, the Dirac monopole singularities are re- 
solved ivitn the fully non-abelian constituent monopoles that appear in the 
calorons,E£3 with_Lh£ only difference that these calorons are now modified by 
the finite volumeE£3 for which no exact analytic solutions are available. They 
nevertheless fit very well to the infinite volume solutions, whereas away from 
the non-abelian cores of the constituent monopoles there is a nearly perfect 
fit to the abelian field of the singular Dirac monopoles. This aspect is illus- 
trated in Fig. 16, comparing the analytic result in the decompactification limit 
with the nuineiical result on a lattice. For further details one should consult 
the literaturellH^I It is natural to conjecture, in the light of these results that 
the fractionally charged instanton is mapped to a pipgle point charge on a 
torus with (abelian) C periodic boundary conditionsEll in the decompactifica- 
tion limit. 



5 Gauge Fields on the three-Sphere 



The most important reason to study gauge fields on the three-sphere is that 
conformal equivalence of S*" 



o IR gives a very simple and explicit con- 
L27| ^Yiig allows one to formulate in a precise 
way how the dependence,, cap be encoded in the boundary conditions on 



struction for the instantons 
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the fundamental domainP^rP^ Due to the curvature of the sphere, at large 



volumes the corrections to the glueball masses are in powers of the inverse 
radius, as opposed to an exponential approach for the torus (to be discussed 
in Sec. 6.1). But ultimately, any geometry can be scaled-up to an infinite vol- 
ume, and should in this limit give the same results. Therefore, by comparing 
different geometries we may indirectly get some useful information. 

We embed 5*^ in by considering the unit sphere (the radius R can be 
reinstated on dimensional grounds where required), parametrized by a unit 
vector n^. Alternative formulations, useful for diagonalizine_lJie Faddeev- 
Popov and fluctuation operators, were developed by CutkoakyJHj We can use 
the 't Hooft tensors r] and fj to define orthonormal framingsli23 of S^, which were 
motivated by the particularly simple form of the instanton vector potentials 
in these framings. The framing for is obtained from the framing of TR'^ by 
restricting in the following equation the four-index a to a three-index j (for 
a — one obtains the normal on S'^, see Eq. (pq)). 



(72) 
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Note that e and e have opposite orientations. Each framing defines a dif- 
ferential operator with associated (mutually commuting) angular momentum 
operators Li and L2: 

It is easily seen that = Lj, which has eigenvalues ^/-|-1), with / = 0, i, 1, • • •. 
By identifying the logarithm of the radius in as timeJn the geometry 
X H, the (anti-)instantons are easily obtained from thosetil on M'' 

_ e^-o-^ _ e.jfccr^'e^ - {v + e^'n^)(J, 

2(l + eMn^)' 2(l + £^n,,) • ^'^^ 

Here e'^, CTj and Aj are defined with respect to the framing for instantons 
and with respect to the framing for anti-instantons, and we introduced 

u=- 7^5 7, = -; r5 s = Aexp(t). (75) 

The unit quaternions cr^ were given in Eq. (p5[). The instanton describes 
tunneling from A = at t = — cxd to = — Cj at t = 00, over a potential 
barrier at t =Jlthat is lowest when 6 = 0. This configuration corresponds 
to a sphaleron^ i.e. the vector potential Aj — —^(Jj is a saddle point of the 
energy functional with one unstable mode, corresponding to the direction (u) 
of tunneling. At t = 00, = —Uj has zero energy and is a gauge copy of 
Aj = by a gauge transformation h — n^a^ with winding number one. 

We will be concentrating our attention to the 18 modes that are degenerate 
in energy to lowest order with the modes that describe tunneling through 
the sphaleron and "anti-sphaleron" . The latter is a gauge copy by a gauge 
transformation h' — n^fj'* with winding number —1 of the sphaleron. The two 
dimensional space containing the tunneling paths through these sphalerons is 
consequently parametrized by u and v through 

A^{u,v) = {~u-el~vel)^. (76) 

The gauge transformation h with winding number one is easily seen to map 
(m, v) = [w, 0) into (m, v) — (0, 2 — w). The 18 dimensional space is defined by 

A,{c, d) = (c,"e-^^ + d,"e^^) ^ = A,{c, d)e^^. (77) 

The c and d modes are mutually orthogonal and satisfy the Coulomb gauge 
condition, djAj{c,d) = 0. This space contains in it the {u,v) plane through 
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c° = ^"'^^ ^^'^ ^"j — — w(5i^The significance of this 18 dimensional space is 
that the energy functional 
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V(c,d) ^ ^ l ltr{F,'^)=V{c)+V{d) + ^{{c1)^d';)^~{c1d'^)'}. 



V(c) = 27T^\2{cfr+6dctc+-Mc-)^-{c-c^)^]\, (78) 



is degenerate to second order in c and d. Indeed, the quadratic fluctuation 
operator M in the Coulomb gauge, defined by 

-^Jtr(F,|) = j^^tv{A,MrjA,)+0{A^), 

- 2Ll5,j+2{Li + S)l, StJ=-iea^J, (79) 

has A{c, d) as its eigenspace for the (lowest) eigenvalue 4. These modes are 
consequently the equivalent of the zero-momentum modes on the torus, with 
the difference that their zero-point frequency does not vanish. 

To find the fundamental region by minimizing the norm functional, Eq. (^, 
we can use FPf(A) in Eq. (|^) as a hermitian operator acting on the vector 
space C of functions h over with values in the space of the quaternions 
m = {g^cr^lq e IR"*}. The gauge group G is contained in C by restricting to 
the unit quaternions: G = {h d C\h = h^^ap^.h e IR'',/i^t/i^ — 1}. This allows 
us to extract detailed information about the Gribov and fundamental region 
within the 18 dimensional space (c, d). When minimizing the norm functional 
over C, instead of G, one obviously should find a smaller space A C A. With 
£ a linear space, A can now also be defined by the condition that FPf{A) be 
positive, 

A = {A e r| {h, FPf{A) h) > 0, V/i e C). (80) 

Similar to the Gribov horizon, the boundary of A is now determined by the 
location where the lowest n on-t rivial eigenvalue of FPf{A) vanishes. For the 
(c, d) space it can be shownEHa that the boundary dk will touch the Gribov 



horizon 957. It establishes the existence of singular points on the boundary of 
the fundamental domain due to the inclusinn A C A C fi. (By showing that 



the fourth order term in Eq. (H) is positivep2 this is seen to correspond to the 
situation as sketched in Fig. g^) 

In order to simplify the notation, write FPi = FPf and FPi = FP, with 
the indices related to the isospin. The associated generators are 

Ti = It and Ti = ladf. (81) 
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We can now make use of the SU(2)xSU(2)xSU(2) symmetry generated by Li, 
L2 and Tt to calculate explicitly the spectrum of FPt{A). One has 

FPMic, d)) = - -^ctT,^L\ - jdtT,^Ll (82) 

which commutes with —1^2, but for arbitrary (c, d) there are in general no 
other commuting operators (except for a charge conjugation symmetry when 
t = i). Restricting to the (u, v) plane one finds that 

FPtiAiu, v)) = 4Li + juU ■ + -^vU ■ T*, (83) 

which in addition commutes with Jt = Li +L2 + Tt, the total angular momen- 
tum, and is easily diagonalized. Fig. ^ summarizes the results for this (m, v) 
plane and also shows the equipotential lines, as well as exhibiting the multiple 
vacua and the sphalerons. As it is easily seen that the two sphalerons are gauge 
copies (by a unit winding number gauge transformation) with equal norm,[tiiey 
lie on 9A, which can be extended by perturbing around these sphaleronsu 
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Figure 17: Location of the sphalerons (large dots), classical vacua (smaller dots), the Gribov 
horizon (fat sections), the boundary of A (dashed parabola) and part of the boundary of the 
fundamental domain (full curves, through the two sphalerons). Also indicated are the lines 
of equal potential in units of 2" times the sphaleron energy. 



To obtain the result for general (c, d) one can use the invariance un- 
der rotations generated by Li and L2 and under constant gauge transfor- 
mations generated by Tj, to bring c and d to a standard form, or express 
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det (^FPt{A{c, d))\i^i^, which determines the locations of dfl and 9A, in terms 

of invariants. We define the matrices X and Y by = c°c^ and = djdj, 
in terms of which 

det (FPi(yl(c,d))|;^i^ = [81 - 18 Tr (X + y) + 24(det c + det d) 

-(Tr(X-r))2 + 2Tr((X-y)2)]^ (84) 

A two-fold multiplicity (the square) is due to charge conjugation symmetry. 
The expression f m f = 1 , that determines the location of the Gribov horizon 



in the (c, d) spacep3 is somewhat more complicated. If we restrict to d = 
the result simplifies considerably. In that case one can bring c to a diagonal 
form = XiS"'. Rotations and gauge transformations reduce to permutations 
of the Xi and simultaneous changes of the sign of two of the Xi . One now easily 
finds the invariant expression ( Tr (X) = J^i ^1 ^-nd det c = Xi) 

det ('fFi(A(c, 0)) 1,^1 ] = (2detc-3Tr(A:) + 27)''. (85) 




Figure 18: The fundamental domain (left) for constant gauge fields on 5'^, with respect to 
the framing e^, in the diagonal representation Aj = XjCTj (no sum over j). The dots on the 
faces indicate the sphalerons. On the right we show the Gribov horizon, which encloses the 
fundamental domain, coinciding with it at the singular boundary points along the edges of 
the tetrahedron. 

In Fig. |l^ we present the results for A and 17. In this particular case, where 
d = 0, A coincides with A, a consequence of the convexity and the fact that 
both the sphalerons (indicated by the dots) and the edges of the tetrahedron 
lie on i9A, the latter also lying on dVl. It is essential that the sphalerons do not 
lie on the Gribov horizon and that the potential energy near 9f2 is relatively 
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high. This is why we can take the boundary identifications near the sphalerons 
into account without having to worry about singular boundary points, as long 
as the energies of the low-lying states will be not much higher than the energy 
of the sphaleron. It allows one to study the glueball spectrum as a function of 
the CP violating angle 9, but more importantly it incorporates for 9 — the 
noticeable influence of the barrier crossings, i.e. of the instantons. 

An effectiveJiamiltonian for the c and d modes is derived from the one-loop 
effective actionE23 To lowest order it is given by 

where g{R) is the running coupling constant (related to the MS running cou- 
pling by a finite renormalization, such that the kinetic term in Eq. ( ^6|) has no 
corrections). The one loop correction to the effective potential is given by 



Viff^(c) = A^i«)2_^«2detc-f K3(c^n' + M(c^?f (87) 



with the following numerical values for the coefficients 

Ki = -0.24534599851796, Kg -0.84996541224534, 

K2 = +3.66869179814223, Kg = -0.06550330854836, 

K3 = +0.50070320309661, K7 = -0.36171221599671, 

Ki = -0.83935963341300, Kg = -2.29535686135471. 

Along the tunneling path (e.g. c° = —uSf and d — 0) the effective potential 
can be calculated with simpler methods, providing an important check. 

Unlike for the torus, where in lowest order all excitations iiLthe zero- 
momentum modes are degenerate and Bloch perturbation theorycfl provides 
a rigorous definition of the effective hamiltonian, one has to rely here on an 
adiabatic approximation not controlled by the coupling constant. The low 
lying excitations in the c and d modes are well below the excitations of t hf 
modes that were integrated out, justifying the adiabatic approximation 
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It provides enough room to achieve a satisfactory understanding of the non- 
perturbative dynamics due to spreading of the wave functional to the boundary 
of the fundamental domain. The boundary conditions are chosen such that 
the gauge and (left and right) rotational invariances are preserved and that 
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they coincide with the appropriate boundary conditions near the sphalerons. 
Projection on the irreducible representations of these symmetries is essential to 



reduce the size of the matrices to be diagonalizec 
All this could be implemented in a tractable way 
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a Rayleigh-Ritz analysis, 
see Fig. 




In perturbation theory the 0~ "oddball" is slightly lighter than the 0+ glue- 
ball state. In lowest order the 0+, 0~ and 2+ glueball states are all degenerate, 
lifted in the next order in perturbation theory, by an amount independent of 
the one loop corrections 

(mo--mo+)i? - -g^{R)/A7r^ + 0{g^), 

(m2+-mo+)i? - 5g2(i?)/67r2 + 0(54), (89) 

whereas m^R = 4 - (1.25 - ki) ^JJl) /2Tr'^ + O (g^). This perturbative re- 
sult was also found by DiekmannE23 (His one loop corrections, however, do 
not survive the test of reproducing the effective potential along the tunneling 
path.) Clearly, the 0~ glueball being the lightest state is an artifact of the 
finite volume, and it is an important test to see if non-perturbative effects 
that set in when going to larger volumes are able to correct for this unwanted 
feature. Indeed, when including the effects of the boundary of the funda- 
mental domain, the 0^/0^ mass ratio rapidly increases from slightly below 1 
to above. At the same time the slow rise of the 2+/0'*" mass ratio becomes 
more rapid. This behavior is observed both with and without the one loop 
corrections included. In the latter case the non-perturbative effects set in at 
g{R) ~ 2. Beyond g{R) ~ 2.8 it can be shown that the wave functional start 
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to feel parts of the boundary of the fundamental domain which the present 
calculation is not representing properly. This value of g{R) corresponds to a 
circumference of roughly 1.3 fm, when setting the scale as for the torus, assum- 
ing the scalar glueball mass in both geometries at this intermediate volume to 
coincide. That instanton effects are largely responsible for the spin splittings is 
further confirmed by studies in the instanton liquid model^^J which concludes 
that instantons lead to an attractive, neutral and repulaiije force in respectively 
the scalar, tensor and pseudo scalar glueball channelsE 
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6 Large Volume Results 

We divide A by the set of all gauge transformations Q, including those that 
are homotopically non-trivial, to get the physical configuration space. All the 
non-trivial topology is then retrieved by the identifications of points on the 
boundary of the fundamental domain. This becomes important when the wave 
functional spreads out in configuration space, which happens at large volumes, 
whereas at very small volumes the wave functional is localized around A = 
and one need not worry about these non-perturbative effects. That these effects 
can be dramatic, even at relatively small volumes (above a tenth of a fermi 
across), was demonstrated for the case of the torus. It has historically been 
very important that the torus with periodic boundary conditions provides a 
relatively wide window, from a tenth of a fermi to three quarters of a fermi, 
where these non-perturbative effects can reliably be taken in to account, and 
were carefully verified in lattice Monte Carlo studies. It shows the intricate 
dynamics due to the intrinsic non-linearities of the theory and the fact that 
the issue of gauge copies can not be ignored and is important for the infrared 
dynamics of the theorju 

The hierarchy of boundary effects is effectively described by the various 
saddle points, that typically lie at the boundary of the fundamental domain. 
An important lesson can be learned from comparing different geometries, be- 
tween which the structure of the fundamental domain deviate considerably. As 
we stressed before, the shape of the fundamental domain is independent of L 
if the gauge field is expressed in units of 1/L. It is only the strength of the 
coupling constant that controls the spreading of the wave functional. Suppose 
that the coupling constant will grow without bound. This would make the po- 
tential irrelevant and makes the wave functional spread out over the whole of 
field space (which could be seen as a strong coupling expansion). If the kinetic 
term would have been trivial, the wave functionals would be "plane waves" 
on a space with complicated boundary conditions. In that case it seems un- 
avoidable that the infinite volume limit would depend on the geometry (like 
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T"^ or S^) that is scaled- up to infinity. Due to the non-triviahty of the kinetic 
term this conclusion cannot be readily made and our present understanding 
only allows comparison in volumes around one cubic fermi. However, one way 
to avoid this undesirable dependence on the geometry is that the vacuum is 
unstable against domain formation. As periodic subdivisions are space filling 
on a torus, this seems to be the preferred geometry to study domain formation. 
It is hard to formulate this in a precise way, let alone find an order parameter 
for this domain formation. But let us assume domain formation does occur. 

Since the ratio of the square root of the string tension to the scalar glueball 
mass shows no structure around L = 0.75 fermi, we assume that within a 
domain both reach their large volume value. The color electric string now arises 
from the fact that flux that enters a domain has to leave it at the opposite 
side. Flux conservation with these building blocks automatically leads to a 
string picture, with a string tension as computed within a single domain and a 
transverse size of the string equal to the average size of a domain. The tensor 
glueball in an intermediate volume is heavily split between the doublet (E^) 
and triplet {T^) representations of the cubic group, with resp. 0.9 and 1.7 
times the scalar glueball mass. This implies that the tensor glueball is at least 
as large as the average size of a domain. Rotational invariance in a domain- 
like vacuum comes about by averaging over all orientations of the domains. 
This is expected to lead to a mass which is the multiplicity weighted average 
of the doublet and triplet, yielding a mass of 1.4 times the scalar glueball 
mass. In the four dimensional euclidean context, 0(4) invariance makes us 
assume that domain formation extends in all four directions. The deconfining 
temperature would again be set by the average domain size. As is implied 
by averaging over orientations, domains will not neatly stack. There will be 
dislocations which most naturally are gauge dislocations. A point-like gauge 
dislocation in four dimensions is an instanton, lines give rise to monopoles and 
surfaces to vortices. In the latter two cases most naturally of the type. 
We estimate the density of these objects to be one per average domain size. 
We thus predict an instanton density of 3.2 fm^**, with an average size of 1/3 
fermi. For monopoles we predict a density of 2.4 fm~'^. The SU(2) spectra in 
volumes around the domain size, for the sphere and the torus, are compatible 
with TO(2+)/m(0"'") ~ 1.5 and m(0^)/m(0+) ^ 1.7. In this picture, we would 
further predict \/K/mo+ ~ 0.24. All these results seem to give the right order 
of magnitude, close to the values one measures in a large volume. 



m(2+)/m(0+) = 1.46 ± 0.09, VK/mo+ = 0.267 ± 0.009, 
m(0-)/m(0+) = 1.78 ± 0.24, (90) 



taken from a recent review by Tepertiia It would be useful to have high- 
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precision lattice Monte Carlo data around L = I fcrmi. In this intermedi- 
ate volume range, the nature of the states changes from being dominated by 
zero-momentum fields to genuine particle states. It is important to take this 
into account in choosing the proper variational basis for the operators in the 
analysis of the Monte Carlo data. It is not sufficiently appreciated that it is 
exactly this change of the nature of the states that can provide fundamental 
physical insight in the formation of the mass gap and the confining string. 

The studies reported so far in small and intermediate volumes were from 
first principles. The volume plays the role of the control parameter to keep the 
strength of the interactions in check. The challenge is to find which are the ef- 
fective degrees of freedom that drive the formation of the mass gap. The lesson 
we learned is that large field fiuctuations are essential. They enter as sphaleron 
configurations in the hamiltonian formulation, and are thus associated with in- 
stantons. Nevertheless, it seems inescapabie-jthat instantons, monopoles and 



vortices all have to be part of the equationP±il in one way or another 



6. 1 Volume Dependence of Stable Particle Masses 

Once a mass gap is being formed, this gives of course the relevant degrees of 
freedom to describe the large distance behavior. For the torus, isolating the 
volume dependence to come from the propagator winding around the bound- 
ary of the boXjphas led Liischer to derive the exponential approach of stable 
particle massescJ The method makes cunning use of the fact that the finite 
volume propagator is given by 1^l{x) = X)n6S«' where A(a::) is the 

propagator in the infinite volume. With a mass gap the sum over periods will 
converge, since A(a;) g^™'^'. Physically A(x + nL) can be identified with 
propagation from to x, going in addition n times around the box. Going 
around only once gives the leading finite volume correction of the order e~*^^. 
The value of M depends on the mass gap, and which type of virtual corrections 
contribute to the self-energy of the particle under consideration. Obviously, 
without interactions their would be no volume dependence. 

Applying these ideas to QCD, one assumes that the long distance behav- 
ior is described by an effective theory. In the presence of light quarks this 
is of course the chiral effective lagrangian. For the calculation of the volume 
dependence of the glueball, it is assumed to be an effective scalar theory. Con- 
finement implies that propagation of a quark around the boundary, unlike in 
a small volume, is not an option for large L due to the confining string it has 
to stretch. Similarly, gluons are confined and thus cannot propagate around 
the volume. It is only the physical colorless states that can propagate around 
the box at large L. Therefore, the fractional periods that appear in the per- 
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turbative expansion with twisted boundary conditions, see Sec. [3.4| , have no 
effect on the finite size corrections in a large volume. Indeed twisted bound- 
ary conditions are devised such that all gauge invariant quantities are strictly 
periodica 

Vertex functions are not affected by the finite volume, and the expression 
of a Feynman diagram in a finite volume is exactly of the same form as in an 
infinite volume, except that integrating over the vertex positions is restricted 
to the finite volume. The exponential suppression is caused by the finite vol- 
ume modification in the propagator, specific for periodic boundary conditions. 
For the sphere, due to the curvature of the manifold, there will be algebraic 
corrections, and the volume corrections are much harder to establish. To keep 
track of the volume dependence universal properties of one-particle irreducible 
vertex functions are studied. A general Feynman diagram I? will correspond 
to the amplitude 

JLiV)=ll n r^dx^{v)exp{zJ2p(.e)-xie))V Y[ AL{x{fii))-xm)), 

(91) 

where p{e) are external momenta. The set of vertices v (wq those that are 
connected to an external line) and propagators € are denoted by V and C The 
product of all vertex factors (coupling constants) is given by V. 

To single out those contributions that correspond to propagating n times 
around the box, one has to take into account the invariance 

x{v) x{v) +m{v)L, n{e)'^n{e) + m{f{e))~m{i{e)), (92) 

where m^{v) e- 77^ and n{£) occurs in the expansion of Al{x). With graph 
theory LiischerE3 showed that the sum over n{£) splits in a sum over the orbits 
[n] and a sum over the integers m{v). Obviously one has dx = dx 

and only the sum over the orbits remains, Jl^J^) = Jl{Pi N), 

^ f^OO 

Jl{VM)= n n/ dx^(i;)exp(z5^p(e)-x(e))V 

X W A{x{f{l)) - x{i{l)) + n{e)L). (93) 
eeC 

The advantage of this decomposition is that for the orbit [n] = [0], this is 
precisely the infinite volume expression, i.e. J'l{'D,[0]) — Jooi'D). When 
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[n] 7^ [0] there is net winding around the boundary, the length of the extra 
winding given by M^([7i]), thereby leading to an exponential suppression, 

JLiV, [n]) « exp[~mLK{{pie)})W{[n])]. (94) 

The constant K{{p{e)}) measures, in units of the mass gap m, the lowest mo- 
mentum flowing through the propagators, constrained by the fixed external mo- 
menta. For a self-energy graph K({p(e)}) = |\/3, with the external momenta 
on-shell p{e) — (im, 0), the only case needed here. Furthermore, = 1 

for a "simple" orbit [1] (where \ui£)\ = 1, for one ling and for all others) and 
W{[n])>V2 for aU other casesB 

To determine the finite volume corrections for the lightest stable one- 
particle state, one takes Lq = oo and Lj — L (this simply implies we always 
have no — 0). The mass is measured using {(j){t)(p{0){L— exp(— M(L)t), where 
(j){t) is an interpolating field for this state. This definition of the mass, also 
used in lattice calculations, coincides with the pole in the finite volume two- 
point function, G^^(iM(L),0) = 0. The two-point function Gl{p) is defined 
by G^^{p) ~ p^ + m? — T,l{p), where the normalization of the self-energy S 

is such that I]oo(j'Ti, 0) = ^^Soo(*"^, 0) = 0, which implies that in leading 

order SM{L) = -^{^l - Soo)(im,0) = -J-r^^){{im,0), [1]). This can be 
expressed in terms of the full four-point Greens function by 

r(2)(p,[l]) = J ^ [6e*«i^G(<z)] G^^\p,q,-p,-q), (95) 

where [6e^'^^^G{q)] comes from summing J2[n]=[i] J l(T^2^ [1]) over the Feyn- 
man diagrams T>2 that contribute to the full two-point function, and making 
use of the cubic invariance, to equate the contributions with rij = ±1 to ni = 1. 

The full four-point function can be related with the help of Schwinger- 
Dyson equations to the one-particle irreducible three- and four-point vertex 
functions. By shifts of integration contours, taking the analyticity domains 
of the vertex functions properly into account, one finds that the dominating 
contribution comes from the on-shell three-point vertex functions, denoted by 
A, a physical three particle coupling constant. With p = {im,0), and both q 
and p + q "on-shell", qi — iy^qj^+Sm^/i and the integral ov£ij q, assuming 
the three point coupling does not vanish, can be shown to giveEj 

SM{L) = - exp(-iV3mL) + O (g-™^) . (96) 

The subleading term, of order e~™^, is due to the contribution coming from 
the four-point function r*^**). 
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The derivation has been given for the volume dependence of the mass gap, 
e.g. for the scalar glueball mass in pure gauge theory. The results also apply 
for the cases of a stable boundstate particle, with mass tjib = 2m — Eb, 
and for a nucleon that couples to a pion. One only needs to adjust the 
"on-shell" condition, qi = i^qj_ + n\ , for the different kinematical situation. 
For the bound state one finds n\ = m? — TOg/4, which leads to 5Mb{L) = 
— 3A^ exp(— /x^L)/(167rm^L), whereas for the nucleon = — m^/4m^, 
giving SMn{L) = — Om^gl^^ exp(— /i^L)/(167rTO^L), where A has been ex- 
pressed in the appropriate dimensionless pion- nucleon coupling constant g^y^. 
For the finite volume correction to the pion mass in QCD, the pion three- 
point vertex vanishes and one gets a slightly more complicated results of pijder 
exp(— TOjri), involving an integral over the forward scattering amplitudeE3 



6.2 Volume Dependence from Scattering Phase Shifts 

Also it was shown ]m' Liischer how to extract scattering phase shifts from vol- 

13q , [l39| 



ume dependenceli^'i!±3 A less rigorous method, based on the notion of pseudc 
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potentials, allows for a transparent way of understanding this all order result 
Consider the reduced hamiltonian H ~ —m^^d^ + V{x.) for two interacting 
particles in the center of mass in ordinary quantum mechanics with a (Bose) 
symmetric reduced wave functions, i^i^x.) — '(/'(— x). When the potential has a 
finite range A {V{x.) = for |x| > A), the wave function is a plane-wave outside 
of the interaction region. Scattering theory tells us there is a unique relation 
between the incoming wave e"*^'^' and the outgoing wave e*'^'^' in terms of the 
scattering phase shift S{k) 

^{x) = e-*'=l"l + e2^«We*'=l"l, \x\ > A. (97) 

These states form a complete basis of scattering states. In a finite volume the 
periodic boundary condition, tplx + L/2) = ip[x — L/2), implies the following 
implicit equation for the momenta 

e^'^C'^e''''^ = 1, (98) 

which holds as long as L > 2A, showing how the volume dependence of two- 
particle states are related to the phase shift. For a vanishing potential the 
phase shift vanishes, S{k) = 0, and one recovers the standard discretization of 
the momenta, k = 2'im/L. 

In more than one dimension a complete set of scattering wave functions is 
given by 

^T2i^) = - tan(<5,(k))n,(|k|r)]r,™(k), (99) 
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where k = k/|k|, r = |x| and ji, ni are the spherical Bessel-functions. The 
energy is given hy E — k.^ /m. It appears that the spherical nature of the scat- 
tered waves makes it impossible to impose periodic boundary conditions. One 
can, however, introduce the notion of a pseudQ4iptential as was introduced for 



the hard-sphere bose gas by Huang and YangElil The essential idea is equally 
simple as powerful. Replace V{'x.) by the simple pseudo potential Vi5(x) such 
that for |x| > A, eq. ((p9|)) is still an exact solution of the relevant Schrodinger 
equation. Subsequently one solves the Schrodinger equation with Vs as its (en- 
ergy dependent) potential, using periodic boundary conditions. We illustrate 
this by making the simplified assumption that all phase shifts vanish, except 
for Sq. In that case one easily verifies that the pseudopotential is given by 

V5(x) = -i!J-tan((5o(k))4(x). (100) 

Since we have to allow for wave functions that are singular at the origin, we 
can extend our class of functions to those that, when averaged over the an- 
gles, have a Laurent expansion X]^-ArC„r". We define j3(x)?/;(x) — co(53(x). 

Alternatively, (53(x)?/;(x) = ^(r-0(x))53(x). We expand the wave function in 
plane waves, suitable for implementing periodic boundary conditions, ip{'x.) = 
L^^ ^ 6„ exp(27rix • n/L). Substituting this in the relevant Schrodinger equa- 
tion ^ ^ 

+ vsi^m^) = — ^(x), (101) 

m ax^ TO 

is easily seen to give 

47rtan(,5o(k)) 

^" = "|k|(k2-(2™A)^) °- ^ ^ 

With Co = X] ^6 thus find a relation between the momentum in the center 
of mass and the scattering phase shift at this momentum 

tan(Jo(k))Zoo(l;q) ^ kL 

27r2|q| ' ^ 27r ^ ^ 

The zeta-function Zoo(s;q) = X^ne^^^'^^ ~ Q'^)^^ is defined through analytic 
continuation in s. 

To ohlain this result cq should be non- vanishing. There can be "singular" 
solutionaHj for which cq = at momenta k — 2Tni/L, if there exists n' e 7L^, 
such that |n| — \n'\. In that case ■(/'(x) = exp(27rm ■ x/L) — exp(27rm' • x/L) 
is regular and vanishes in the origin. If ■0 is in the scalar representation (Ai) 
of the cubic group, this singular behavior only occurs when n' is not related 
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to n by a cubic transformation, which will make the momentum where this 
occurs quite large. Furthermore, restricting to the Ai sector, eq. (( |103[ )) will 
also be valid if the phase shifts only vanish for angular momenta ^ > 3. This is 
because a spin 2 wave function decomposes in the E and T2 representations of 
the cubic group and hence does not couple to the scalar sector (note that due 
to the Bose symmetry all odd angular momentum phase shifts will vanish). 

In field theory it can be shown that the reduced hamiltonian is replaced by 
an effective Schrodinger equation, that can be derived from the Bethe-Salpeter 



equation for the four-point functionE£3 Its energy-dependent "potential" is pro- 



portional to the so-called Bethe-Salpeter kernel, with the range A determined 
by the polarization cloud (which is why we restrict the analysis to field theo- 
ries with a mass gap) . The fully relativistic two-particle energy W is given by 
W = 2\/k2"+"r7? = 2^m(ra + E). where k is the center of mass momentum 

all order analysis is based on studying the 
- k^)?/'(x) = in a finite volume, 



of the scattering pair. Liischer'gHB 
solutions of the Helmholtz equation (9^/9x^ 
allowing for power-like singularities at the origin. He demonstrates that trun- 
cating to a finite number of phase shifts Sg, {£ < (-max) in general converges 
rapidly with Imax- For the simplified case that only Sq is non- vanishing, the 



result is as in Eq. (103). 

One of the applications is computing the energy of the two-particle states 
as a function of L. Using 



Zoo(l;q) = - 



1 



Zoo(l;0) + q2Zoo(2;0) + 0(g4) 



tan((5o(k)) = ao|k|+C'(fc^ 



(104) 



introducing Zoo{s: q) = Zoo(s; q) — (— q^) and oo, t he s o-called the scattering 
length. One can now straightforwardly iterate eq. ((103)) and find 



Cl 



Ana, 



L3 (1 

^oo(i;0) 



L 

C2 = 



L2 



^oo(i;o) 



^00 (2; 0) 



(105) 



The numerical value£B'Elll are Ci = -2.837297 and C2 = 6.375183. When the 
two-particle energies come close to the mass of a resonance, like the p ov K 
meson coupling to two-particle pion states, the^hase shifts get large and it 
becomes imperative to use the all order resultsE 
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Many of the aspects of the large volume expansions for the one-particle 
and two-particle masses-Jjiave been verified— sjiccessfuUy for the two dimen- 
sional non-linear sigmalll3 and Ising modeltl2l, and for the four dimensional 
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phaseElJ Recently the method 



for studying on the lattice non- 



Isin^dS and 0(4) model in the symmet 
has been revisited by Lellouch and Liischcit 
perturbatively the decay K tttt, important for understanding direct CP 
violation. By studying the decay in a finite volume, the transition amplitude 
of the kaon to the discrete set of finite volume two-pion final states can be re- 
lated to the infinite volume decay rate. With the advance in computer power 
this calculation will be feasible in the future. 



6.3 Goldstone Modes and Chiral Perturbation Theory 

The finite size corrections to the string tension K{L) also rely on an effec- 
tive description, namely that of the bosonic relativistic string. The finite size 
behavior is described in terms of the force by 

with K the infinite volume string tension. The so-called Liischer termEli with 
the 1 / L"^ correction is universal. A simple way to understand this correction is 
in terms of the Casimir energy in one space dimension (the string) on a periodic 
interval of length L (the length of the string) for two massless scalar fields 
(the two independent transverse fluctuations of the string), 2 • i 2'n:\k\/ L — 
47rC(— 1)/L (in dimensional regularization, with C, the Riemann zeta function). 
Liischer's original derivation was for a string with fixed end-points, for which 
the Casimir energy is a quarter of the one with periodic boundary conditions. 
The transverse fiuctuations are examples of Goldstone bosons, that here arise 
due to the spontaneous breakdown of the transverse translational invariance 
in the background of a stringlilZ 

Goldstone bosons do of course also play a very important role in the low- 
energy description of QCD in the form of pions. They occur with massless 
quarks due to the spontaneous breakdown of chiral symmetry. With an explicit 



mass term for the quarks, the pions acquire a mass and the results of Sees. 3.1 



and 3.2 apply as long as m^L ^ 1. For massless pions this can of course not 
be realized, and instead chiral perturbation theory becomes an effective tool. 
It also applies for ^ provided L is in the so-called mesoscopic range, 
which is rii^^ 3> i S> Ag^^. The principle behind chiral perturbation theory 
is that the symmetries of the low-energy theory strongly constrain the effective 
lagrangian. The resulting sigma model has its symmetry determined by the 
number of light quark flavors and the gauge group. For SU(3) with two light 
flavors, the up and the down quark, this is the 0(4) sigma model spontaneously 
broken down to 0(3) by the chiral condensate. Finite size effects occur in 
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powers of 1/L, but are determined by infinite volume quantities^omputed 



essentially through replacing the continuous by discrete momentacl^ 



From the practical point of view this gives access to these infinite vol- 
ume quantities through finite size effects. An important example is the chira. 



condensate, which can be related through the Leutwyler-Smilga sum ruleslll^ 
to eigenvalue distributions of the Dirac operator (averaging over the gauge 
fields). It is the average spacing of these eigenvalues near zero that through 
the Banks-Casher formulaE£2l gives access to the chiral condensate (AA) = 
7r(L'^T|(i/)?/))|)~^. The roughness of the gauge field causes (AA) to be inversely 
proportional to the volume (L^T) of space-time. At weak coupling, like in a 
small volume, the typical eigenvalue spacing would be but in addition 

we saw in Sec. B.2 there are no near zero eigenvalues and chiral symmetry is 
manifest. 

Remarkably, in the mesoscopic domain the low-lying eigenvalue distribu- 
tion {pLXX) is related to a universal function called the microscopic spectral 
densityl^ Psiu) = limv^ooiVE)^^ {p (^{VT.)^^u^) . It can be calculated in 
Random Matrix Theory, leaving the chiral condensate_E = \{'4'4')\ cis a single 
free parameter. We refer to a review by Verbaarschott£3 for details. 



6.4 Electric- Magnetic Duality 

Much of the finite volume work in gauge theories has started with 't Hooft's 
paper on the torus with twisted boundary conditionsB We come back to it one 
more time to show how it leads to an electric-magnetic duality, which can be 
used to derive finite size results for the energy of magnetic fiux in terms of the 
string tension, assuming one has a confining string. 

Consider the gauge group SU(A^) and fix the magnetic fiux m by twisted 
boundary conditions. In the ^0 = gauge this can be specified as in Eq. (pl|). 
Given these boundary conditions, the remaining gauge freedom is specified by 
the twisted gauge transformations /ik, defined in Eq. (^5[), where k classifies 
the element of the homotopy group H2(T^ , Zn) ~ Zf^. Two twisted gauge 
transformations with the same value of k S Zf^ differ by a gauge transformation 
with integer winding number. A particular representative of can be chosen 
such thatE2l vihk) = — k-m/A^, as we argued below Eq. (pO|). The generalization 
of Eq. ( |2C| ) to m ^ can now be written as 

^{[h^h^,] A) = e^''^ exp (^^^^^ ^(A), e, ^ e - |-m. (107) 

The non-trivial dependence implies that magnetic fiux will carry in general 
a fractional electric flux. This result is intimately related to the observation of 
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Wittenty that the 9 dependence can be implemented at the hamiltonian level, 
by replacing the canonical momentum belonging to the gauge field, i.e. the elec- 
tric field -E'i(x), by i?i(x) — -^Bi, with Bi{x) = ^eijkFjk{x) the magnetic field. 
This behavior, of magnetic charge obtaining a 9 dependent electric chaige, also 



plays a role in 't Hooft's abelian projection and oblique confinementE£i 
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Figure 20: The SU(6) spectral fiowCj of Eg in a small volume, Eq. (108), representing the 
energy of electric flux for m = (0, 0, 1) and e\ = 
period 2it, the flow only with a period of 12tt. 



0. The spectrum is periodic with a 



It seems that the 2tt periodicity in 9 needs to be replaced by 27r7V. Nev- 
ertheless, the spectrum itself is periodic with period 2tt. There is, however, 
a non-trivial spectral flow, which in the infinite volume may be non-analytic 



and give rise to oblique confinement The spectral flow can neatly be if 



lustrated in a finite volume, where energies are required t o b e analytic in 
9. For this we consider m = (0, 0, 1) and recall from Sec. 3^ that the en- 
ergy splittings for e = (0, 0, 63) are due to tunneling with fractionally charged 
instantons. One ea sily 



Stanton calculations 



imputes the resulting energies through standard in- 
using the fact that the twisted instanton has only 
four zero-modes due to translations (cmp. Sec. m), 
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-S^^) =Csm^{Tr\eg\/N)exp[-87r^/Ng\L)]/Lg\L). (108) 



The behavior of Eg is illustrated, for = 6, in Fig. from which it is 
immediately clear that the spectrum is periodic with a period of 27r, but that 
following the energy levels adiabatically, one only returns to the same situation 
after N = 6 periods. 
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We now consider the euclidean partition function at finite temperature, 
specified by f//3, wtiere /3 is tfie period in the euclidean time direction. For 
simplicity we will also assume the period lattice S is orthogonal, a''^-' = L^e^^^'^ 
with Lq = (3. The free energy F{e, m, L^) is defined for the different superse- 
lection sectors as 

exp [-(3Fg{e, m, L^)] - Tr,^^ [exp(-/3i7^)] , (109) 

where the trace is here over the Hilbert space with definite 9, electric and mag- 
netic quantum numbers. To relate the finite temperature partition function to 
the euclidean path integral, 9 and the electric flux have to be implemented as 
a sum over the homotopy classes (y ^TL and k e Z^) of the gauge transfor- 
mations. In the Aq = gauge, the gauge field at i = is related to that at 
t = /3 by ^(x, /3) — [/ik/ii] ^(x, 0). For a given choice of v and k we denote 
the euclidean path integral by M^,y(k, m, L^) — W^(n, S), where n denotes the 
antisymmetric integer twist tensor, with riy — tijimi and rtoi = fci specifying 
the twisted boundary conditions, which can now be formulated in any gauge. 
The relation between the free energy and the euclidean path integral reads 

exp[-/3Fe(e,m,L^)] =iV-3 ^ ^ e-«z;''-«iy,(k, m, L^), (110) 

where we introduced the A''-th root of unity, z^^ = e^^^/^ . 

The alternate notation, Wj/(n^i^, S), makes manifest that from the point of 
view of the euclidean path integral there is no essential distinction between the 
twisted boundary conditions in the space and the time directions. As observed 
by 't Hooft, this path integral is invariant under a joint S0(4) rotation of 
the period lattice S and the twist tensor n. This nearly trivial fact, however, 
results in a duality between the electric and magnetic sectors.EI A simultaneous 
rotation in the 1-2 and the 3-0 planes over 90 degrees, interchanges (fci, fe) with 
(mi,m2), leaves and ma u nchanged and maps tOiL^ = {L3, L2, Li, /3). 
Substituting this in Eq. (|llO| ) leads to the exact dualityH 

exp[-/3Fe(e,m,L^)] = N-^^S^,^^Jg^^^ z|™~™ '= exp[-L3Fe(e, m, Z^)]. 

(111) 

It is because eg • m — ra • eg, is independent of 9, that 9 does not explicitly 
appear in this equationllj Hence, we will in the following ignore the 9 label. 

Let us normalize F(0,0, L^) = 0, and call those fiux states for which 
Lfj,F{e, m, L^) tends to zero light fluxes. That they cannot all be light is easily 



69 



seen by summing Eq. (Ill) over mi, m2, ei and 62 



N-^ exp[-/3i^(e,m,L^)] =exphL3i^(e3e(3),m3e(3),i^)]. (112) 



If (eae^'^^ TOae^'^)) produces a light flux state, there must be iV^ — 1 additional 
light flux states (with the same 63 and 7713). On the other hand, since the 
euclidean path integral is positive, irrespectiv e the twist, once the flux (e, m) 
is light, so must be the flux (0,m) (see Eq. (|lTo|)). For SU(2) and SU(3) it 
can thus be concludeda that, in the conflned phase, the magnetic fluxes have 
a vanishing free energy in the inflnite volume limit. 

It will be assumed that for f3 ^ Li the free energy factorizes in a magnetic 
and electric component, F{e,m., L^) = Fm(m, L^) + Fe{e,Li^i). This can be 
justified by the fact that the electric flux is squeezed into a string with a flnite 
width (to give rise to a linear potential), which at high (3 takes negligible 
space, whereas magnetic flux spreads out over the whole volume (for the free 
energy to vanish). We now show how to compute the volume dependence of 
the magnetic energy, £',„(m, L^) = lim^^oo(-F'm(ni, L^) - F,„(0,L^)), for the 
case where 7713 — 0. Let us choose 63 = and (for simplicity) Li = L2 = L, 



with /?, L3 > L. From Eq. (Ill) it follows that 



Cexp[-/3F™(mi,L^)] =iV-2 J2 z^^'™^ exp[-i3-Fe(gi, i^)]- (113) 

The normalization constant C (obtained by putting m = 0) can be absorbed in 
Fm ■ One can compute the electric free energies from the statistical distribution 
of flux strings of energy KL, running in either the 1- or 2-direction, to build up 
the fluxes from rii+Ci strings running in the i-direction and in the opposite 
direction. For SU(2), fluxes running in opposite directions are equivalent, and 
one instead sums over rii even for Ci — and odd for = 1. Entropy 
allows us to ignore higher units of electric flux (including the ones running 
diagonally). Strings running in the long 3-directions can be neglected as well. 
The Boltzmann weight of one string equals LL3 exp{—/3KL)/A, where is a 
constant related to the "proper area" of the string. It is now straightforward ti 



compute the electric free energy, substitute the result in Eq. (113) and obtai 



the magnetic energy by extracting the limiting behavior in /3 00, 

E„, 2{2-SN:2)[sin^iTrmi/N) + sin^{TTm2/N)]A~^Lexp{-LL3K). (114) 

Thus the magnetic energy falls off in leading order with the same area law as 
the Wilson loop expectation value. 
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7 Conclusions 



I have tried to convince the reader of the usefulness of a finite volume as a 
control in studying four dimensional non-abelian gauge theories, with for the 
torus the practical benefit of providing a guide for lattice Monte Carlo results. 
There are many issues I have not addressed here. Needless to say the choices 
have been determined by the things I have worked on in the past, or that were 
close to me. I hope to have at least succeeded in developing this theme here 
in a logical and pedagogical fashion. It is my belief that probing the non- 
perturbative dynamics through spectral properties is the only reliable way to 



gain insight. This is not to say it is only through finite vo 
make progress. There is of course the large A'^ expansion 
for a string representation of QCD 
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tudies we may 
and the hope 
It should also be said we have followed 
the standard paradigm, that confinement is first to be understood in the pure 
gauge sector, instead of taking_the attitude that the light quarks provided by 
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nature is what mostly matters 

On the basis of this review we list a number of open problems that are 
worthwhile addressing. 

1. Solve the problem with the adiabatic approximation in the calculation 
of the Witten index in small volume supersymmetric gauge theories. 

2. Find analytic solutions for the basic instantons on the torus, with and 
without twist. 

3. Isolate the associated sphaleron configurations and non-perturbatively 
important degrees of freedom. 

4. Map out with high-precision the low-lying spectrum around 1 fcrmi. 

5. Establish in a reliable way the physical nature of the finite volume cross- 
over near 1 fermi, and find the intermediate-distance effective degrees of 
freedom, that are responsible for the mass gap and string formation. 

The first two items are interesting technical problems in their own, which one 
may hope are close to being solved. They have a wider range of applica- 
tions than in the context of QCD alone. The fourth item does not have to 
wait for analytic predictions, and can be addressed by Monte Carlo methods 
with present day resources. The last item is of course what it is all about. 
It is unlikely this can be solved by a single mathematical equation, or with 



mathematical rigorE£j Fortunately, we do have experiment to guide us here. 
Despite the many phenomenological models, victory can not be claimed yet. 
The scenario of Mandelstam and 't Hooft, for describing QCD in terms of a 
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dual superconductorP3li£3 remains an appealing one and gets some pjjpport 
from supersymmetric gauge theories exhibiting Seiberg-Witten dualityp£j even 
though their matter content is uncomfortably far from that of QCD. 

The challenge of really understanding QCD non-perturbatively remains 
wide open, even though many in the particle physics community escape to 
other dimensions. Some of them hope that from this higher vantage point we 
can look down upon this problem that was left over from the twentieth century, 
and solve it with twenty-first century techniques. Anything is allowed, as long 
as we come up with a practical and reliable solution. 
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have all contributed in an essential way to developing the material presented. 
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